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COEFFICIENTS 



WOLTER GROENEVELT 



Abstract. The irreducible ^-representations of the Lie algebra su(l, 1) consist of discrete se- 
ries representations, principal unitary series and complementary series. We calculate Racah 
coefficients for tensor product representations that consist of at least two discrete series repre- 
sentations. We use the explicit expressions for the Clebsch-Gordan coefficients as hypergeomet- 
ric functions to find explicit expressions for the Racah coefficients. The Racah coefficients are 
Wilson polynomials and Wilson functions. This leads to natural interpretations of the Wilson 
function transforms. As an application several sum and integral identities are obtained involving 
Wilson polynomials and Wilson functions. We also compute Racah coefficients for U q (su(l, 1)), 
which turn out to be Askey- Wilson functions and Askey- Wilson polynomials. 



1. Introduction 

In this paper we study Racah coefficients, or 6j'-symbols, for representations of the Lie alge- 
bra su(l, 1). The Racah coefficients for certain tensor product representations of su(l, 1) lead 
to unitary integral transforms with a very-well-poised 7-Fg-function, called a Wilson function, 
as a kernel. These Wilson functions and the corresponding integral transforms were recently 
introduced by the author in |12j with the applications in this paper in mind. We also consider a 
g-analogue of su(l, 1), namely the quantized universal enveloping algebra U q (su(l, 1)). From the 
Racah coefficients for certain tensor product representations of U q (su(l, 1)) we obtain a new in- 
terpretation of the Askey- Wilson functions transform introduced by Koelink and Stokman \2'A\ . 
As an application we obtain several identities for special functions involving (Askey-) Wilson 
functions and polynomials. 

Racah coefficients for su(2) play an important role in the theory of angular momentum in 
quantum physics p]. They were first studied in the 1940's by Racah |32j . who also obtained an 
explicit expression as a finite single sum for these coefficients. Only much later it was realized 
that the Racah coefficients are multiples of polynomials of hypergeometric type, so that the 
orthogonality relations for the Racah coefficients lead to discrete orthogonality relations for 
certain polynomials. These polynomials are nowadays called the Racah polynomials, which can 
be defined by 

, , . f-n,n + a + f3 + l,-x,x + j + 6 + 1 

R n (x = R n (x; a, f3, 7, 6 = 4^3 ,1 a , x . 1 11 ;1 

\ a + 1,(3 + d + 1,7 + 1 

where one of the lower parameters is equal to —N, N € Z>o, and < n < N. These are 
polynomials of degree n in the variable x{x + 7 + 5 + 1), and they are orthogonal on the set 
{0,1,..., A}; 

JV 

^2 Rm(x)R n {x)W(x) =0, m^n, 

x=0 
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where 



W(x) = W(x;a,p,-y,5) 



2 7 + 25 + 2x + 1 (a + l) x Qg + £ + l) g (7 + l) x ( 7 + 5 + 1) : 
2 7 + 25 + l ( 7 + 5 - a + l) x ( 7 - /? + l) x (5 + l) x x! 



From their definition it is clear that the Racah polynomials satisfy the duality property 



so the Racah polynomials are also polynomials of degree x in n(n + a + /3 + 1). This shows 
that the dual orthogonality relation for the Racah polynomials are essentially the same as the 
orthogonality relations. In section [2] it is shown how the orthogonality relations for the Racah 
polynomials can be obtained directly from their interpretation as Racah coefficients. We refer 
to the lecture notes by Van der Jeugt jHH] for a detailed treatement of how properties of Racah 
polynomials can be obtained from their interpretation as Racah coefficients, see also Vilenkin 
and Klimyk |39| . 

The Wilson polynomials, introduced by Wilson jlU], can be defined by 



Assuming that a,b,c,d > these polynomials satisfy continuous orthogonality relations on 
[0, oo), with the notation f{a ± b) = f(a + b)f(a — b), 



so we may consider the Wilson polynomials as continuous extensions of the Racah polynomials. 

Together with the Racah polynomials the Wilson polynomials form the top level of the Askey- 
scheme of hyper geometric polynomials |2U| . This is a large scheme of orthogonal polynomials in 
one variable that can be written as hypergeometric series. All the polynomials below the top 
level in the Askey-scheme, and many of their properties, can (at least formally) be obtained by 
taking appropriate limits of the Wilson or Racah polynomials. Therefore these polynomials are 
fundamental objects in the theory of hypergeometric orthogonal polynomials. 

In the Wilson functions were introduced. The Wilson function can be considered as 
an analytic continuation of the Wilson polynomial in its degree, or equivalently as continuous 
extension of the Racah polynomials in both the variable x and the dual variable (the degree) 
n. The Wilson funtions satisfy generalized orthogonality relations, in the sense that they are 
the kernel in two different unitary integral transforms, called the Wilson transforms I and II. 
The dual integral transforms are essentially the same as the integral transforms, which can be 
considered as continuous analogues of the fact that the orthogonality relations and the dual 
relations for the Racah polynomials are essentially the same. 

The main goal of this paper is to show that the Wilson functions have a natural interpretation 
as Racah coefficients for the Lie algebra su(l,l). Both Wilson function transforms appear in 
this context, and unitarity of both integral transforms can be obtained in this way. Moreover, 
we show that the Wilson polynomials also appear as Racah coefficients for su(l,l). As an 
application we obtain several sum or integral identities for Wilson polynomials and Wilson 
functions. 

Let us also mention some other generalizations of the Racah polynomials. Askey and Wilson 
j5] considered basic hypergeometric extensions of the Racah polynomials, and called these the q- 
Racah polynomials. The g-Racah polynomials satisfy orthogonality relations on a finite discrete 
set, similar to the Racah polynomials. Later Kirillov and Reshetikhin gave an interpretation 
of the g-Racah polynomials as Racah coefficients for the quantized universal enveloping algebra 
Z/^(su(2)). Askey and Wilson 4_ showed that the g-Racah polynomials also satisfy continuous 
orthogonality relations, and considered as a polynomial in a continuous variable, these poly- 
nomials are the now famous Askey- Wilson polynomials. Koelink and Stokman [23] introduced 



R n {x; a, P, 7, 5) = R x (n; 7, <5, a, f5) 
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the Askey- Wilson functions, which can be considered as the analytic continuation of the Askey- 
Wilson polynomials in their degree, and they studied the corresponding integral transform. In 
section E] we show that the Askey- Wilson polynomials and the Askey- Wilson funtions have an 
interpretation as Racah coefficients for U q (su(l, 1)). 

Let us now describe the contents of this paper. In section El we recall the irreducible *- 
representations of the Lie algebra su(l, 1). The irreducible representations consist of four classes: 
positive and negative discrete series tt^, principle unitary series ir p and complementary series 
tt c . Then we calculate Clebsch-Gordan coefficients and Racah coefficients related to tensor 
product representations of positive discrete series. The results from section |2] are well known 
in the literature. This section is mainly included to introduce notations and to explain the 
ideas we use later on. Also we use the results from section |2] in following sections. In section 
we first decompose 2-fold tensor product representations, of which one is a discrete series, into 
irreducible representations, and we calculate the corresponding Clebsch-Gordan coefficients. 
The results from this section are used later on in sections 01 El El to compute Racah coefficients 
related to various 3-fold tensor product representation of su(l, 1). In sectionEJwe compute the 
Racah coefficients for the tensor product 7r + (g> ir~ tt + . This leads to a natural interpretation 
of the Wilson transform I. In section E] we consider Racah coefficients for the tensor product 
7r + (gi 7r p (g) 7T~, and in this way we obtain a natural interpretation of the Wilson transform II. 
In section El we give the Racah coefficients for some other tensor products involving at least two 
discrete series representation. In section [7| we use natural identities for Racah coefficients to 
obtain sum and integral identities involving Wilson polynomials and Wilson funtions. Finally 
in section El we sketch how the Askey- Wilson polynomials and the Askey- Wilson functions can 
be obtained as Racah coefficients for U g (su(l, 1)). 

We use the standard notations for (basic) hypergeometric functions as in the book by Gasper 
and Rahman jllj . 

Acknowledgement. I would like to thank Erik Koelink for many useful discussions and sug- 
gestions. 



2. REPRESENTATIONS OF su(l, 1) AND RACAH POLYNOMIALS 

This section is meant to explain the ideas we use later on in the paper to find Racah coeffi- 
cients, and to introduce some notations. 

2.1. The Lie algebra su(l,l). The Lie algebra su(l, 1) is generated by H, E and F, which 
satisfy the commutation relations 

[H, E] = 2E, [H, F) = -2F, [E, F] = H. 

There is a *-structure defined by H* = H and E* = —F. The center of the universal enveloping 
algebra U(su(l, 1)) is generated by the Casimir element fi, which is given by 

Q = ~^(H 2 + 2H + 4FE) = ~^(H 2 - 2H + 4EF). (2.1) 

The irreducible ^representations of su(l,l) are determined by Bargmann [3], see e.g. also |391 
§6.4]. There are, besides the trivial representation, four classes of irreducible ^representations 
of 5u(l,l): 
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Positive discrete series. The positive discrete series representations ir^ are labelled by k > 0. 
The representation space is £ 2 (Z + ) with orthonormal basis {e n } ng z> - The action is given by 

n£(H)e n = 2(k + n)e n , 



tt+(E) e n = ^{n + l)(2fc + n) e n+1 , 

ir£(F)e n = - y/n(2k + n - I) e„_i, 

7r A |"(^)e n = fe(l - k)e n . 

Negative discrete series. The negative discrete series representations ir^ are labelled by k > 0. 
The representation space is ^ 2 (Z>o) with orthonormal basis {e n } n6 z> - The action is given by 

n^(H)e n = -2(k + n)e n , 



7r k (E)e n = - y / n(2k + n - 1) e n _i, 

n' k (F)e n = V(n + l)(2/c + n) e n+1 , 

vr^(Jl) e n = fe(l - fe) e n . 

Principal unitary series. The principal unitary series representations 7r^ e are labelled by 
e G [0, 1) and p > 0, where (p, e) 7^ (0, |). The representation space is ^ 2 (Z) with orthonormal 
basis {e n } n£ z- The action is given by 

TTp (H)e n = 2(e + n)e n , 



Kp, £ ( E ) e n = J{n + e + ^-ip)(n + e + ^ + ip) e n+1 , 



j> - 1 • 1 



(2.4) 



n p ,e( F ) e n = -\/{n + £- --ip)(n + e--+ip)en-i, 
rf, £ (ty e« = (p 2 + 7) e„. 



For (p,e) = (0, ^) the representation ir^ splits into a direct sum of a positive and a negative 
discrete series representation: 7r^i = 7r 
subspaces: {e n j n < 0} {e n | n > 0}. 



discrete series representation: 7r.f1 =7Ti©7t 1 . The representation space splits into two invariant 

' "2 ~ ~ 



Complementary series. The complementary series representations 7r£ e are labelled by e and 
A, where e G [0, |) and A G (—55 — e ) or e G (|, 1) and A G (— ^,e — 1). The representation space 
is £ 2 (Z) with orthonormal basis {e n } ne z- The actions of the generators H, E, F are the same 
as in the principal unitary series, with ip replaced by A + \. 

We remark that in case e g" [0, 1) the principle unitary series iTp £ are unitarily equivalent 
to the -7Tp £+p , where p G Z is such that e +p G [0,1). This follows from the fact that for 
e G + 1), p G Z, and for X G su(l, 1), the action of tt^ £ (X) on the basis {e n } ng ^ is the 
same as the action of ^ e _ v acting on the basis {e n + p } rag z- A similar observation can be made 
for the complementary series. 

The operators defined above are unbounded closable operators, with domain the finite linear 
combinations of the basis vectors. The representations are ^-representations in the sense of 
Schmiidgen 31. Ch.8]. 

2.2. Clebsch-Gordan coefficients for ir + ® tt + . We consider the tensor product representa- 
tion 7r^ (g) 7rjjT . Recall that for the definition of the a tensor product representation we need the 
coproduct. For a Lie algebra element X the coproduct A is defined by A(X) = 1 <S> X + X ® 1. 
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The coproduct can be extended to the universal enveloping algebra as an algebra morphism, 
and then we obtain from (|2.1|) for the Casimir element 

A(fi) = l<g)0 + Oigil- -H®H-F®E-E®F. (2.5) 

In order to decompose the tensor product 7r^ ® 71^ into irreducible representation we diago- 
nalize the operator 7rjjT ® 71"^ (A(f2)). Applying this operator to e n ®e p _ n G ^ 2 (Z>o)® 2 and using 
(|2.2|) and q2.5|) , we find an expression of the form 

^fei ® 7r fe 2 (^(^)) e ™ ® e p-« = °™ e «+i ® e p-«-i + K e n tg> e p _ n + a n _ie„_i ® e p _ n+ i, 

for some coefficients a n and 6 n . From writing out a n explicitly we see that a_i = a p+ i = 0. So 
for p G Z>o the space 

ft p = span c {e n ® e p _ n | n = 0, . . . ,p} ^ C p+1 

is invariant under the action of the Casimir operator. Let v p (X) = X^ = o v n(A)e n ® e p-n G W p 
be an eigenvector of tt^ ® 7r^(A(0)) for eigenvalue A, then the coefficients t>n are solutions of 

A<(A) = a n < +1 (A) + 6 n <(A) + a^t^A). 

This shows that the coefficients Vn can be generated from Vq. To give an explicit expression for 
the coefficients Vn we introduce the dual Hahn polynomials. 
The dual Hahn polynomials are defined by 

T n (x; 1 ,5,N) = 3^(" n '"^ 1 + 2^ + 1 ;l), n = 0,l,...,N, 

where N is a non-negative integer. These are polynomials in x(x + 7 + 5 + 1). Note that the 
dual Hahn polynomial can be obtained as a limit case of a Racah polynomial; T n (x; 7, 5, iV) = 
lirng^oo R n (x; —N — 1, /3, 7, 5). The dual Hahn polynomial is not symmetric in 7 and 5, but we 
can interchange 7 and 5 using 

T n (x; 7 ,<5,iV) = (-ir^^T N . n (x;5, 7 ,N), (2.6) 

(7 + l)x 

where we assume that x G {0, 1, . . . , N}. The orthonormal dual Hahn polynomials are defined 

by 

t n (x)=t n (x; 7 ,5,N) = {-ly ^^^ ^T^ 7, S, N), 
and they satisfy, for 7 > — 1 and 5 > — 1 or for 7 < —N and 5 < —N, 

N 

^2t m (x)t n (x)W(x;-f,5, N) 

x=0 

= (2x + 1 + 5 + 1)( 1 + 1U-NU5 + 1)n 

{ ,7 ' ' 1 (-ir(x+ 7 +s+i) N+1 (s+i) x x\ ■ 

The dual Hahn polynomials satisfy the following three-term recurrence relation, 
-x(x + 7 + 5 + 1) t n (x) = a n t n+1 (x) + b n t n (x) + a n _it n _i(x), 

where 

a n = yj{n + l)(n + 7 + l)(n - 5 - N){n - N), 
b n = (n -I- 7 + l)(n - JV) +n(ra - <5 - JV - 1). 

Note that for x G {0, ...,N} the polynomial T n {x) is also a polynomial in n of degree x. 
Considered as a polynomial in n, T n {x) is called a Hahn polynomial. 
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We see that the coefficients Vn are multiples of dual Hahn polynomials with parameters 
(7,(5, N) = (2ki — 1, 2&2 — l,p), and the eigenvalue A is given by A (x) = (fei+fe 2 +x)(l — k\ — k 2 — x). 
We normalize the coefficients Vn as follows. We split the weight function for the dual Hahn 
polynomials in a p-dependent and a ^-independent part, such that the p-dependent part equals 
1 for p = 0, and we glue the square root of the p-dependent part to the dual Hahn polynomial. 
So for p G Z>o we set 



vlix-MM = \ \S\^1 X( 1 2) \ P t»(x;2fci-l,2fc2-l,p), x = 0,l,...,p. 
y {2ki + 2k 2 + x) p 

Let VF(x; fei, ^2) be the weight function given by 

w( u u \ {2x + 2k l + 2k 2 -l)(2k 1 ) x 

W(x; kx,k 2 ) = t — — w r — r, x G Z> . 

(x + 2fci + 2fe 2 - l)(2A;2)a;ic! 

We denote by C p+1 (VF(x; fei, ^2)) the (p + l)-dimensional Hilbert space with the inner product 
associated to this weight function at the points x G {0, 1, . . . ,p}. Then, for fixed p, the functions 
x h- ► Un(x; fei, /C2) form an orthonormal basis for C p+1 (iy(x; fei, ^2))- We remark that it is 
customary to multiply the coefficients Vn by the square root of the weight function W, since 
then they are orthogonal in C p+1 with respect to the usual inner product. We will work instead 
with weighted inner products. We can now formulate our finding as follows. 

Proposition 2.1. For p € Z>o, the operator I : 7i p — ► C p+1 (W(x; fei, k 2 )) defined by 

I : e n <g) e p _„ m|ih u£(x; fei, k 2 j\ , 
is unitary and intertwines n k <S> ft^ (A(U)) with M^ kl+k2+x ^i_ kl _ k2 _ x y 

Here M a is denotes multiplication by a. We have found that 

p 

V ' P { X ) = v n( x 'i h, k 2 ) e n ® e w _ p 

n=0 

is an eigenvector of ir'fc <g) ^ k2 ( A (Q) ) for eigenvalue ( fei + fe 2 + x ) ( 1 — fei — fe 2 — x ) . The squared norm 
of i> p (x) is W{x;k\,k 2 ) . We write u p (x) = v p (x; k\, k 2 ) if we need to stress the dependence 
on the representation labels fei and fc2- A comparison of the multiplication in Proposition 12.11 
with the action of f2 in the discrete series representation (|2.2|) and (|2.3j) . suggests that for all 
x G Z>o the discrete series representation 7r kl+k2+x appears in the decomposition of ir~£ <8> 7r^ . 
To find the exact decomposition we determine the action of the generators H, E and F on the 
eigenvector v p (x). The action of H is easy to find, and the actions of E and F follow from the 
following contiguous relations: 

(d - 1) 3*2 (;-_\ C e ; l) = a 3 F 2 (° + ^ C ; l) + (« . - i + 1) 3 F 2 ; 1 

W-W-c),*^.! 1 )- W 

rf(e - a) 3^2 ( d ^ 5 1 J + d(a + b + c - d - e) 3 F 2 I ^ ^ : 

The first relation is obtained from combining (3.7.9) and (3.7.11) from [2], the second relation 
from combining (3.7.10), (3.7.12) and (3.7.14) from |2j. This gives us the following actions 

7T+ <g> 7T+ (A(#)) U P (X) = 2(fel + fe 2 + P ) V P (X), 

ir+®ir+ 2 (A{E))v p (x) = y / (p-x + l)(2fei + 2fe 2 + x + p) v p+1 {x), 
7T+ ®ir+(A(F))v p (x) = y/{p - x)(2fei + 2fe 2 + x + p - 1) ^(x). 
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Observe that these are the same as the actions of the generators in the representation ^ +k2+x 
on the basis vectors e p - x . Therefore we define 

oo 

I' : £ 2 (Z> )® 2 -> ®l 2 (Z>o)W(x;k l ,k 2 ), 

x=0 
ni+712 

e ni ® t n2 I > ^ v nl^ 2 ^1 ' ^2) e rii+7i2— xi 
x=0 

which corresponds to inverting v p (x) = ^ n v p l (x; ki, k 2 ) e n (8 e p - n using the orthogonality rela- 
tions for the dual Hahn polynomials. Here the weighted direct sum can be considered as a direct 
integral with a measure \x supported on Z>q, satisfying ^({x}) = W(x; k%, k 2 ) for x £ Z>o- Now 
we have the following proposition. 

00 

Proposition 2.2. The operator I[ intertwines 7r^ 18 ir^ with ^^^x+ks+x' 

Using the notation tti = H2 if the two representations Hi, i = 1,2, are unitarily equivalent, 
Proposition 12.21 states that the tensor product representation 71^ <8 7rjT is decomposed into 
irredubible representations as 

00 

j=0 

and the corresponding Clebsch-Gordan coefficients are multiples of (dual) Hahn polynomials. 
See also |25| . |27j . |35j for the interpretation of (dual) Hahn polynomials as Clebsch-Gordan 
coefficients. 

2.3. Racah coefficients for ir + (8 7r + <8 7t + . To determine the Racah coefficients for the tensor 
product n ki <8 7r^ 2 (8 7rjf we first determine the eigenvectors of Q,. We consider the operator 
n ki ® ^fca ® ^fe ® Using the actions of H, E, F on v p (x), we obtain that 7r^ (8 71^ (g> 

7r^" 3 (A <8> 1)(A(0)) acts on v p ~ m (x) (8 e m in the same way as ir kl+k2+x 8> vr^" 3 acts on e p _ x _ m (8 e m . 
So we see that the vector 

p— XI 

W P ~ X1 (xo) 



p—Xl p—U3 

= ^2 ^ Vn~ n3 (xi;k 1 ,k 2 )v P l ~ Xl (x2;k 3 ,k 1 + k 2 + x\) e m ® e n2 <8e re3 , 



n,3=0 ni=0 

where ni + ^2 + 713 = p, is an eigenvector of ir k (8 7r k2 (8 55 l)(A(f2)) for the eigenvalue 

(k\ + /c2 + &3 + xi + x 2 ){l — k% — k 2 — k% — x\ — x 2 ). Moreover, the actions of H, E, F, on Xl (x 2 ) 
are the same as the actions of the generators in the representation TT kl+k , 2+k . i+Xl+X2 on e p - Xl - X2 . 

The squared norm of the vector w I [~ Xl (x 2 ) is [W{x\] ki, k 2 )W(x 2 ; k 3 , k\ + k 2 + x±, )) 1 . 
In the same way we find that the vector 
p-yi 

w p f vi {y 2 ) = < yi (y 2 ; h,k 2 + k 3 + yi ) e ni ® ^- ni (yi; h, k 2 ) 

ni=0 

P-Vi p-ni 

= ^ X] v n 3 " ni (yi; fc 3,^2X7 !/1 (y2;A;i,fc2 + £3 + yi)e ni <8e„ 2 (8 e„ 3 , 

ni=0 ri3=0 

where p = n\ + n 2 + n 3 , is an eigenvector of 7r^ (8 vr^ (8 7r^ (1 (8 A)(A(Q)) for the eigenvalue 
(ki + k 2 + k 3 + y 1 + y 2 )(l - ki - k 2 - k 3 - y x - y 2 ). 
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From the associativity of the tensor product, or equivalently the coassociativity of the co- 
product, i.e. (A <g) 1)A = (1 <g> A)A, it follows that we have found two bases of eigenvectors for 
the same operator. Let us consider the overlap coefficients for the transition from one basis to 
the other, i.e. the coefficients 

C PljP2 (xi,x 2 ; yi ,y 2 ;h,k 2 ,k 3 ) = {w p l 1 ' Xl {x 2 ),w p 2 2 ~ yi (y 2 )). (2.8) 

Here the (•, •) denotes the usual inner product in ^ 2 (Z>o) <x ' 3 . Since both w\ and w 2 are eigenvec- 
tors of ft, the coefficient C vanishes if x\ + x 2 ^ y\ + y 2 - Also, from applying H to w\ and using 
H* = H, we see that the coefficient vanishes if p\ ^ p 2 . Similarly, setting pi = p 2 = p, applying 
E on w\ and using E* = —F, it follows that C is independent of p. From here on we denote, 
Cp, p {xi,x 2 ; y±,y 2 ; h,k 2 , k 3 ) = U(x\,yi\ k\, k 2 , k 3 ,t), where we assume that t = x\ + x 2 = yi + y 2 . 
The coefficient U is called the Racah coefficient. 

Let us define a weight function in the variable x as follows; 

W(x; fci, k 2 ,k 3 ,t) = W(x;ki, k 2 )W(t - x; k 3 , k x + k 2 + x, ). 

So this weightfunction is the inverse of the squared norm of w p x {t — x). Note that interchanging 
k\ and k 3 gives the inverse of the squared norm of w p ~ x (t — x). Now we have 

t 

w p 2 ~ y {t -y) = U (x, y; h,k 2 , k 3 ,t)w p ~ x (t - x)W(x; k\, k 2 ,k 3 ,t), 

x=0 

t 

w[~ x (t-x) = ^C/(x,y;A;i,A;2,A;3,t)'u;f" s/ (t - y)W(y;k 3 ,k 2 ,k 1 ,t). 

This gives us the following orthogonality relations for the Racah coefficients 

t ^ 
/l U (x, y; h,k 2 , k 3 ,t)U(x, y r ;kx,k 2 , k 3 ,t)W(x; h,k 2 , k 3 ,t) - 

x=0 
t 



(2. 



W{y;k 3 ,k 2 M,ty 



y^U(x,y;k 1 ,k 2 ,k 3 ,t)U(x',y;k 1 ,k 2 ,k 3 ,t)W(y;k 3 ,k 2 ,k 1 ,t) = — 
^ W(x:k 1 ,k 2 ,ks,t) 

Also note that there is an obvious symmetry property; 

U(x,y;ki,k 2 ,k 3 ,t) = U{y,x\k 3 ,k 2 ,ki,t). 

In order to find an explicit expression for the Racah coefficient, we take the inner product 
with e ni ® e„ 2 <g> e n3 , ni + n 2 A- n 3 = p, in the first identity in l|2.9|) . This gives 

v p n - ni (y; h, k 2 )v p ; y (t -y;k 1 ,k 2 + k 3 + y ) = 

_ ( 2 -10) 

} y £/'(a;,y;fci,fc2,/^ l t)^ 1 ™ 3 (x; ki,k 2 )v p 3 x (t - x;k 3 , ki + k 2 + x)W(x; ki, k 2 , k 3 ,t). 

x=0 

Using the orthogonality relation for the Clebsch-Gordan coefficient v p ^ n3 this formula gives an 
expansion of the Racah coefficient in Clebsch-Gordan coefficients. Now an explicit expression 
for the Racah coefficient can be found by specializing 723 and p in a suitable way. This leads to 
expression as a terminating balanced 4i*3-series, and then we find that the Racah coefficient is 
a multiple of a Racah polynomial, 



TT( 1. fr u *\-r i Y x+ y -t V-(2k 2 ) x (2k 2 ) y I (2fci + 2k 2 + 2fc 3 + t) x {2h + 2k 2 + 2fc 3 + t) y 
U[x,y,ki,k 2 ,K 3 ,t) —{ Ij {2kl ) x{ 2k 3 ) y )j (2k, + 2k 2 + x) x (2k 2 + 2k 3 + y) y 

x R y (x; 2k 2 - I, 2k 3 -1,-t- 1, 2k\ + 2k 2 + t - I). 
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See [38] for details. The above orthogonality relations for the Racah coefficients lead to orthog- 
onality relations for the Racah polynomials on the finite set {0,1,... , iV}. 

In the following sections we study other tensor products for su(l, 1) and the corresponding 
Racah coefficients. In the cases we consider we cannot give eigenvectors of the Casimir operator, 
but instead we must work with generalized eigenvectors. Therefore the computations will be 
more complicated, but underlying ideas are much the same as in this section. 



3. CLEBSCH-GORDAN COEFFICIENTS 

In this section we decompose various 2-fold tensor product representations involving discrete 
series representations, and we give explicit expressions for the corresponding Clebsch-Gordan 
coefficients. These decompositions and the corresponding Clebsch-Gordan coefficients are well- 
known in the literature, see e.g. [18] . In this section we emphasize the relation between the 
Clebsch-Gordan coefficients and orthogonal polynomials. 



3.1. Clebsch-Gordan coefficients for tt + ®it . First we consider the tensor product tt + ® tt 
acting on £ 2 (Z> )® 2 . We define for fixed pGZ elements ft € £ 2 (Z> )® 2 by 

[e n+p <g>e n , p > 0. 

So {/n | n € Z>o, p 6 Z} is an orthonormal basis for ^ 2 (Z>o)® 2 . We define a subspace 
U\ C £ 2 (Z> f 2 by 



U\ = span c {/£ | n G Z> } r(Z> ). 

From ()2.5|) . (|2.2j) and (|2.3j) we find that the operator irt <S> tt^ (A(O)) acting on elements f„ 
extends to an unbounded Jacobi operator acting on TL\. We refer to the book by Akhiezer [Q for 
more information on Jacobi operators. The Jacobi operator we have here is the Jacobi operator 
corresponding to the three-term recurrence relation of the continuous dual Hahn polynomials, 
cf. Prop.2.1]. 

The continuous dual Hahn polynomials are polynomials in x 2 defined by 

/ —n,a + ix,a — ix \ 
b n {x;a,b,c) = (a + b) n {a + c) n3 F 2 [ ;1 . (3.2) 

\ a + o, a + c J 

The continuous dual Hahn polynomials can be obtained as a limit of the Wilson polynomials 
[4Uj . see also [2] or [2U]. By Kummer's transformation [5] Cor. 3. 3. 5] the polynomials S n are 
symmetric in a, b and c. For real parameters a, b and c, with a + b, a + c, b + c positive, or for 
one real parameter and a pair of complex conjugates with a positive real part, the continuous 
dual Hahn polynomials are orthogonal with respect to a positive measure. The orthonormal 
continuous dual Hahn polynomials are given by 

, s / , s (-l) n S n (x;a,b,c) 
s n (x) = s n (x;a,b,c) - 



^n\(a + b) n (a + c) n (b + c) n 



Assuming without loss of generality that a is the smallest of the real parameters, the polynomials 
s n are orthonormal with respect to the measure fJ>(-;a, b, c), given by 

/-y roo _^ 
f(x)dfi(x;a,b,c) = — f(x)w(x;a,b,c)dx+ ^ f (i(a + k))w k (a;b,c), 



a+k<0 
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where 



w(x; a, b, c) 
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r(a ± ix)T(b ± ix)T(c =L ix) 



r(a + 6)r(a + c)r(6 + c)r(±2ix) ' 

r(b - a)r(c - a) (-l) fc (2a) fc (a + l) fc (g + b) k (a + c) fe 
r(-2a)r(6 + c) fc!(a) fc (a-&+l) fc (a-c+l) fc 



Wk(a;b,c) =i Res w(x;a,b,c) 

x=i(a+k) 



The polynomials s n form an orthonormal basis for the Hilbert space C{p) consisting of even 
functions (//-almost everywhere) that have finite norm with respect to the inner product associ- 
ated to the measure p. The orthonormal continuous dual Hahn polynomials satisfy the following 
three-term recurrence relation, 



x s n (x) (x) + b n s n (x) + a n _is n _i(x), (3.3) 



where 



a n = \/{n + l)(n + a + b)(n + a + c)(n + b + c), 

6 ra = 2n 2 + 2n(a + 6 + c — -) + ab + ac + be. 

The Jacobi operator we obtained from the action of the Casimir corresponds to the three-term 
recurrence relation for the continuous dual Hahn polynomials with parameters 



{a,b, c} 



{ki - k 2 + §, ki + k 2 - |, fc 2 - ki - p + i}, p < 0, 
{/c 2 - fci + |, fei + fc 2 - |, fci - k 2 +p + \}, p > 0. 



The corresponding orthogonality measure depends also on p. Using F(a +p) = (a) p T(a) for the 
continuous part of the measure, and (a + p)k = (a)k( a + k) p /(a) p for the discrete part, we split 
this orthogonality measure in a p-dependent part and a p-independent part, and we glue the 
square root of the p-dependent part to the continuous dual Hahn polynomials, i.e. for p € [0, 00) 
we define the functions F% by 



F%{p-M,k 2 ) 



(k 2 - h + \ ± ip)- p 



V (-p)!(2fc 2 )_ p 
xs n (p; ki - k 2 + \,kx + k 2 - \,k 2 - ki - p + \), p < 0, 



(-1) 



n /(A;i-A: 2 + 5±«p)p 



V p! (2A;i)p 

X«n(p; &2 - &i + 5, fci + fc 2 - |, fci - fc 2 +p + i), p > 0, 



and the functions are defined similarly for p in the discrete part of the support of p\, where 
pi denotes the corresponding orthonormality measure. Explicitly, if k\ and k 2 are such that the 
measure is absolutely continuous, p\{-;ki,k 2 ) is given by 



wi(p; h,k 2 ) 



/f'OO 
f(p)dpi(p;ki,k 2 ) = f(p)w 1 (p;k 1 ,k 2 )dp 

1 r(fci - k 2 + \ ± ip)T{k! + k 2 - \ ± ip)T(k 2 -k l + \±ip) 



2vr 



T{2k 1 )T{2k 2 )T{±2ip) 



for all p € Z. We note that in general the number of discrete mass points of p± may still depend 
on the value of p, and that this number is finite for all p. Since the weights of the measures p\ 
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do not depend on p, we can define the measure p\ to be the measure with support 

(J supp(/^(-;fci,fc2)), 



is unitary and intertwines 7r^ (g> 7r fc2 (A(il)) wrai/i M 



and with the same weight function as all the other measures p\. 

Remark 3.1. In order to make the situation as simple as possible we assume without much loss 
of generality that k\ and k 2 are such that the measure fj^ does not have discrete mass points. 
In this case the measures p\ are the same for all p, so p\{-\ k±, k 2 ) = pi(-; fei, fo). 

Proposition 3.2. For p € Z i/ie operator I 2 : — ► £(pi(p; ki, k 2 )) defined by 

-p 2 +\- 

Observe that for p > 0, p 2 + | corresponds to the action of the Casimir in the principal 
unitary series. 

Remark 3.3. Applying the operator I\ to e n ® e m gives F™~ m or F^~ m , depending on the sign 
of n — m. To simplify notations, we use both expressions regardless of the sign of n — m. This 
is justified by the following identity for continuous dual Hahn polynomials, cf. |14[ (3.13)], 

Sn(p;h - k 2 + -, h + k 2 - -, k 2 - h - p + -) = 

(— I) 53 (fci - k 2 + - ±ip)pS n - p (p;k 2 - k\ + -,h + k 2 - -,h - k 2 + p - -), p> 0. 
From Proposition 13.21 we find that 

oo 

v[{p) = Y,n{p-MM)z, 

n=0 

is a generalized eigenvector of 7r^ ® 7r^" 2 (A(r2)) for eigenvalue p 2 + |. We can also determine the 
actions of H,E,F on v^(p), but, since v\{p) £ 2 (Z>o), this will only give formal expressions. 

Let {e n } ng z denote the standard orthonormal basis of £ 2 (Z). We define the unitary operator 
/(by 



/^ 2 (Z> f 2 ^ £ 2 (Z)dp 1 (p;k 1 , i 



o 

"OO 



/■oo 

fn^ F%(p;k 1 ,k 2 )e p dp 
Jo 



Then, under the assumption that pj^ does not have discrete mass points, we have the following 
proposition. 

oo 
r® 

Proposition 3.4. The operator I[ intertwines 7r^ (g) tt^ 2 luzf/i / ^hi—ki^P- 

o 

Here we use that, for e € [p,p+ 1) for some p € Z, the representation 7rf" E is equivalent to the 
representation 7T^ e _ p . 

To proof the proposition we only need to check that /{ indeed intertwines the actions of the 
generators H, E and F. The action of H follows from a straightforward computation. The 
actions of E and F follow from expressing the continuous dual Hahn polynomials as 3i*2-series, 
and using the contiguous relations (|2.7|1 . We refer to (141 Thm.2.2] for a detailed proof. 
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In case the measure pi(-; ki, k 2 ) has discrete mass points, these discrete terms correspond to 
a finite number of discrete series representations, or exactly one complementary series, cf. |141 
Thm.2.2]. 

3.2. Clebsch-Gordan coefficients for ir + ®7t p . Next we decompose the tensor product 7r + ® 
ir p and we compute the corresponding Clebsch-Gordan coefficients. This can be done in the 
same way as the tensor product 7r + ® 7r~, i.e. by diagonalizing the action of the Casimir. 
We define for fixed p G Z the vector c$ G ^ 2 (Z> ) <g> ^ 2 (Z) by 

5n = e n ® e p _ n , (3.4) 

and we define the subspace by 

= span c {^ I n G Z> } = ^ 2 (Z> ). 

From (|2.5[) . I|2.2|) and (|2.4[1 it follows that the action of the Casimir element in the tensor 
product 7r^!~ (g> ir p £ acting on elements extends to an unbounded Jacobi operator acting on 
Ti.2- This Jacobi operator corresponds to the recurrence relation for continuous dual Hahn 
polynomials with parameters 

{a, b,c} = {k + ip,k - ip,- - k - p - e}. 

As before, we split the orthonormality measure for these polynomials in a p-independent part 
and a p-dependent part, and we glue the square root of the p-dependent part to the continuous 
dual Hahn polynomial, i.e. for r G [0, 00), we define 



G p n (r;k, P ,e) 



"I" g — ~tt" s n {r; k + ip,k - ip,\ - k - p - e), r G [0, 00), 

J + k + e ± ZTjp 2 



i + £ ± ip)„ 1 

s n {Tj;k + ip,k-ip,--k-p-E), r = Tj e iR <0 . 



I (j + l) p (2k + 2e - j) p 

(3-5) 

Here we denote Tj = i(i — k — e + j), where j G Z is such that Tj G iM<o- Observe that it follows 
from the explicit expression as a 3i*2-function, that p 1— > Gn(T;k,p,e) is a continuous function 
for /) G [0,oo). For fixed p G Z, the functions G« form an orthonormal basis for the Hilbert 
space C(p^(r\ k, p, e)j , where is the measure defined by 

/roc 
f(T)di4(T-,k,p,e) = f(T)w 2 (T]k,p,£)dT+ V f(T j )w 2 (T j ;k,p,e), 
J ° <irz7 

Jfc£.>_p 

Tj<0 

where 

1 T(k + ip± ir)T(k — ip± ir)T{\ - k - e±ir) 

w 2 (t: k, p, e) = — - — — , rG 0, 00), 

V H ' ' 2tt T(2k)T{\-e±ip)T(±2i T ) V 1 

_ T{2k + e-\±ip) (-i)j(i- 2 fc- 2e) j $-k-e) j $-e±ip) j 
w 2[ r, K ,p,e) T(2k + 2e _ l)T(2k) j\(l-2k-e±ip) j {\-k-e) j ' T Ty 

Remark 3.5. The number of discrete mass points of pF 2 tends to infinity as p — > 00. 

Let p 2 (-; k, p, e) be the measure correspond to the above defined weight functions w 2 and with 
support 

IJ supp(p?,(-;k,p,e)). 



So, by Remark 13.51 this measure has an infinite number of discrete mass points. 
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Proposition 3.6. For p € Z the operator I 2 : — * ^(/ U 2( T ' ^ P> e )) defined by 
is unitary and intertwines tt^ ® 7r^ £ (A(Q)) wif/j M r2+ i 



So we find that 



n=0 

is a generalized eigenvector of 7r^ (g> 7r^ e (A(fJ)) for eigenvalue t 2 + \. We define the unitary 
operator by 



4 :£ 2 (Z> )®£ 2 (Z) -> J £ 2 {Z)w(T;k,p,e)dT®^e 2 {Z> )w(T j ;k,p,e), 

J 

/oo 
GP(T;k,p,e)e p dT + ^ G p n {Tj]k, 
Pi e ) e p-j ■ 

3&>- P 

TJ<0 

Then we have the following proposition. 

00 
/•© 

Proposition 3.7. The operator I' 2 intertwines ® £ with I ir^/dr © ^t'+j' w ^ iere 

J 

e' = e + and f/ie direct sum is over all j £ Z suc/i i/iai e' + j > 0. 

The proposition is proved by checking that intertwines the actions of the generators, using 
the contiguous relations (|2.7|) . The explicit expression for the Clebsch-Gordan coefficients as 3-F2- 
series can also be found in ^H], but the continuous dual Hahn polynomials are not mentioned 
there. 

We remark that the tensor product -k^ <g> tt^ £ can be treated in completely the same way 
as irt <8> T^oe- Formally, the Clebsch-Gordan coefficients can be obtained from the substitution 



k "p,e 

p h-> — i(X + |). 

3.3. Clebsch-Gordan coefficients for tt p (g) ir~ . Let 1? : su(l, 1) — > su(l, 1) be the involution 
defined on the generators of su(l, 1) by 

= -H, ${E) = F, = E. 

This involution can be extended to Z^(su(l, 1)) as an algebra morphism. The Casimir element 
satisfies #(Q) = Q. If we identify the representation space for with the representation space 
for 7T^~ , we have nt ($(X)) = 7rr(Jf) for X G U(su(l, 1)). Furthermore, define a unitary operator 
U : £ 2 (Z) -» £ 2 (Z) by C/e n = (-l) n e_ n . Using J23J it is easy to check that £/o7r£_ £ otf = vr^of/, 
so 7r^_ e o -d is equivalent to 7r^ £ . 
For p E Z, we define 

Z£ = e n+P ®e n G^(Z)®^ 2 (Z>o), (3.6) 
or equivalently = (— l) n +P Ue- p - n ® e n . Compare this last expression with (|3.4|) . From the 
above considerations we obtain 

< £ 0vr-(A(X))C = (-ir +p <_ £ 0vr+(AWX)))5- p , X G W(fiu(l, 1)), 

where g% is the vector g% with the factors in the tensor product interchanged. So we can obtain 
the decomposition of tt p £ ®t:^ from the decomposition of 7rjjT <£> ftp- £ (recall that the coproduct A 
is cocommutative). Since = f2 and since the Clebsch-Gordan coefficients are eigenfunctions 
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of f2 in the tensor product, we obtain from Proposition 13.61 that the Clebsch-Gordan coefficients 
for the tensor product 7r^ £ <g) 7r^ are the functions ( — l) n+p Gn P {&] k, p, — e). From the tensor 
product decomposition for 7r^~ (g) 7r^ e we obtain 

oo 
r@ 

rf,e ®^k~ ^e,da TTj_ e , , £ ' = £ - k, 

i 
where the direct sum is over all j £ Z such that j — e' > 0. 

Finally let us mention that the tensor product of two principal unitary series can be de- 
composed as a direct integral over the principal unitary series, plus an infinite sum of discrete 
series. However, the principal unitary series occur with multiplicity 2. For this reason we avoid 
tensor products of two principal unitary series in this paper. See for instance JJj for related 
Clebsch-Gordan coefficients. 



4. RACAH COEFFICIENTS FOR 7T+ (g> 7T (g) 7T+ 

In this section we find explicit expresssions for the Racah coefficients related to the tensor 
product 7r + (gi tt~ 7r + , and we relate the corresponding generalized orthogonality relations to 
the Wilson function transform of type I. 

4.1. Racah coefficients. We consider the tensor product representation 7r^ (g> 7r^ (g> 7r^, and 
we want to find the corresponding Racah coefficients. First we consider the action of the Casimir 
Q, in the tensor product (n^ <8> vr^) <8> w^, and we want to diagonalize this action. We find that 
the subspace of ^ 2 (Z>o)® 3 spanned by elements e m <8>e n2 <g>e n3 for which n\ — n 2 + n 3 is constant, 
is invariant under the action of the Casimir. Therefore we define, for r € Z, 

where fn is defined by Then, for r 6 Z, the subspace 



X = span c {A£ m |n,m € Z> } C ^ 2 (Z> )® 3 , 

is invariant under the action of £1. 
We define 

K n,m( T ,P;h,k 2 ,k 3 ) = F^~ m (p;ki,k 2 )G r m (T;k 3 ,p,k 1 - k 2 ). 
Let rrf be the product measure defined by 

f(T,p)dm r (T,p;k 1 ,k 2 ,k 3 ) = j j f(T,p)d(p r 1 ~ m (p;k 1 ,k 2 ) x /i£(r; k 3 , p, h - k 2 )). 

The set \K r nm \ n ^ m ^L >a is an orthonormal basis for H(m r (x,y;ki,k 2 ,k 3 )^, the Hilbert space 
consisting of functions in two variables, even in both variables, with inner product associated to 
the measure m r . 

Proposition 4.1. For r £ Z the operator J\ : A r — > 7^(m r (r, fcj, /c2 5 ^3)) defined by 

J l : fc n,m -» ((^P) ^ K, m ( T iP; k li k 2,k 3 )^J, 

is unitary and intertwines 7r^ <g> 71"^ (g> 7r^((A (g> 1) o A(f2)) M r2+ i . 

Proof. First observe that Ji is unitary since it maps an orthonormal basis onto another or- 
thonormal basis. 
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For the second statement we first apply the operator I[ from Proposition 13.41 in the first and 
second factor of the tensor product: 

I[ ® id : k r n m = f n - m ® e m ^ / F^- m (p; ki, k 2 )e r - m ® e m dp. 



o 



Here id denotes the identity operator on £ 2 (Z>o). From the intertwining property of I 2 , see 
Proposition 13. 4| we find 



I 2 ® id o 7T+ vr fc2 7T+ ((A ® 1) o A(fi)) = / tt+ (A(O)) dp. 



Next we apply I 2 from Proposition 13.61 then we find that 

/ffi / /-oc \ 

Ji <ip : k r n m i— > ( r i— > / F^~ m (p; k 1} k 2 )G T m (T; fe 3 , p, fe a - fc 2 ) dp j 

o ^ ° ' 

intertwines 7r^ ® ^ ® 7r^((A ® 1) o A(f2)) with M r2+ i. Clearly M r2+ i is independent of p, 

so the operator J\ also intertwines the action of f2 with M T 2 , i, for almost every p. Since the 

function K r nm {r,p) is continuous in p, the proposition follows. □ 

From Proposition 14. II it follows that 

oo oo 

w\{p,t) = G r m (T;k 3 ,p,ki - k 2 )v^' r {p) ® e m = ^ K n,m( T iP'> k U k 2, h) k r n m 

m=0 n,m=0 

is a generalized eigenvector of 7r^ ® ir^ ® 7r^ g ((A ® 1) o A (17)) for eigenvalue r 2 + |. 

Combining Propositions 13 . 2l and 13 . 6l we can construct an explicit intertwiner J{ to decompose 
n ki ® ^fe ® * n *° i rre ducible representations. So J{ intertwines 7r^ ® 7rj~ ® 7r^ with 



oo / oo \ 

/ / < £ '^©0^ +i )dp, 

V J / 



where e 1 = k\ — k 2 + k 3 . We abuse notation and write the intertwiner as 

J'\ ■ Ki,m^> jj K n,m( T iP; k li k 2,k 3 )e p dTdp, 

So e p inside the integral is a basis vector of £ 2 (Z) for r € [0, oo) and e p should be replaced by 
ej-p € £ 2 (Z>o) when r is in the discrete part of the support of the measure dm. 

Next we consider the action of the Casimir in tensor product 7r^ ® (tt^ ® vr^), and we want 
to have a unitary operator J 2 similar to J\ in Proposition 14.11 Recall that the coproduct A is 
cocommutative, so we can use Proposition 13.41 to decompose 7r^ ® 7r^. Then we can find J 2 in 
the same way as above by considering the action of f2 on the element k r n m € A r , where we write 

k n,m as 

k r - p 6d f r ~ n 

and fn is the same as fn defined by (|3.1j) . but with the factors interchanged, i.e. 

e n - p ® e n , p< 0, 



tv 

J n 



e n ® e n+p , p>0. 



By the symmetry in k\ and k 3 for the tensor product tt^ ® 7r fc2 ® 7r^ , we can find J 2 directly 
from Proposition 14.11 bv interchanging k\ and k 3 , and n and m. 
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Proposition 4.2. Let r £ Z i/ie operator J 2 : *4 r — ► 7i{m r (T, a; k 3 ,k2, defined by 

J 2 ■ K,m (i T i cr ) ^ Kn.n^^^fa^fa) 
is unitary and intertwines tt^ 7^ ^((1 A) o A(Q)) u>ii/i Af r2+ i . 
So 

00 

n,m 

is a generalized eigenvector of 7r^ (8) vr^ 7r^" ((1 A) o A(J1)) for eigenvalue r 2 + |. Again we 
can make an intertwining operotar J' 2 similarly to J[. 

From the coassociativity of the coproduct it follows that there exists a unitary operator that 
intertwines i rt 7rj~ o (1 A) o A with 7r^ 71^ 71^ o (A 1) o A. Then it follows from 
Propositions \£l\ and H31 that there exists a unitary operator 

U : H(m r (T,p;fa,k 2 ,k 3 )) -> TC(m r (r' , a; k 3 , k 2 , fa)) , 

such that U o J\ = J 2 . We write £7 as an integral operator, 

(Uf)(r',cr)= // f{T,p)U r T , {T,p;fa,k 2 ,k 3 )dm r {T,p;fa,k 2 ,k 3 ), 



for some kernel C/L (r, p; fci, k 2 , k 3 ). We call this kernel the Racah coefficient for the tensor 
product 7T^ 7r^~ 7r^. We justify this terminology in Eemark 14,41 We will determine an 
explicit expression for the Racah coefficient, using 



K m,n( T 'i a 'i k 3,k2,fa) = j j K^ m (r, p]fa,k 2 ,k 3 )U^ p]fa,k 2 ,k 3 ) dm r (t, p;fa,k 2 ,k 3 ), (4.1) 

which follows directly from U o J\ = J 2 . First we determine some useful properties of the Racah 
coefficient. 

Since both J\ and J 2 intertwine the action of Vt with M 2.1, U leaves functions in the variable 
t fixed, and this gives 

V r T iAr,p;fa,k 2 ,k 3 ) = 0, if r / r'. 

Let v r {p; fa, k 2 , k 3 , r) be the measure m r (t, p; fa, k 2 , k 3 ) restricted to a line on which r = con- 
stant, i.e. 

f(p) dv(p;fa,k 2 ,k 3 ,T) = j f(p)w 1 (p,fa,k 2 )w 2 (T;k 3 ,p,fa - k 2 ) dp. 

For fixed r this measure du r does not depend on r, therefore we will omit the superscript r. 
Then, instead of using (|4.1jl . we can determine an explicit expression for the Racah coefficient 
from 

/>oo 

K m,n( T ^ cr i k s,k2, fa) = / K' r nm {T,p-,fa,k 2 ,k 3 )U r TCT {T,p-,fa,k 2 ,k 3 )dv(y;fa,k 2 ,k 3 ,T). (4.2) 

JO 

We use this identity for the following lemma. 

Lemma 4.3. The Racah coefficient does not depend on r. 

Proof. The functions K r mn (r, a; k 3 ,k 2 ,fa) are the "matrix elements" of the unitary operator J' 2 
that intertwines 7T^ 7r^ tt^ with J* J® 7r^ e drda, where e = fa — k 2 + k 3 . So from the action 
of E we find 

AK m,n-i( T , o-;k 3 ,k 2 ,fa) + B K T m ^ n {r, a; k 3 ,k 2 ,fa) + C K r mn (r, a; k 3 ,k 2 ,fa) 

1 ■ , irr+l, , , , x ( 4 ' 3 ) 

= |r + e + - +zr|ii^ n (T,cr;A;3,K2,A;i), 
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where we assume for simplicity that r € [0, oo) and we denoted 



A = yjn(2ki +n), B = yjm(2k 3 + m- 1), 
C = \/ (m + n — r)(2k2 + m + n — r — 1). 

Applying (|4.2|) this gives, 



\r + e+ - + ir\ I U^ 1 (r,p;k 1 ,k 2 ,k 3 )K^J l {T,p;k 1 ,k 2 ,k 3 )d^(p;k 1 ,k 2 ,k 3 ,T) 
z Jo 

= J (AK^_ 1)m (r, p;h,k 2 ,k 3 ) + B K T n ^ m _ x (r, p;k 1 ,k 2 ,k 3 ) + C K r n ^ m (T, p\ h,k 2 , k 3 

x U^ a (r, p;k 1 ,k 2 ,k 3 )di^(p;k 1 ,k 2 ,k 3 ,T) 
1 f°° 

= \r + e + - + ir\ / Ul a (T,p-k lj k 2 ,k 3 )K^(T,p-k ll k 2 ,k 3 )du(p-k 1 ,k 2 ,k 3 ,T). 
1 Jo 

In the last step we applied (|4,3|) again with (n, k\) <-> (m, /C3). From the above identity it follows 
that the Racah coefficient is independent of r. □ 

Remark 4.4. Let us explain why we call the kernel U^, (r, p; ki, k 2 , k 3 ) the Racah coefficient. 
Formally, the Racah coefficient can be defined in terms of the generalized eigenvectors by 

U' (p, a; k l7 k 2 , k 3 ,T)5 rir2 5 T y = {w r 1 1 (p,T),w r 2 2 {a,T'))p {z ^ Q)(S i. 

Compare this expression with (|2.8|) . So the Racah coefficients can be considered as "matrix 
elements" of the operator U' defined by U' : w\{p^t) h-> w 2 (a,T 2 ), and as before the coefficient 
is independent of r. If we "expand" one eigenvector in terms of the other basis of eigenvectors, 
we have, for p 2 = y, 

w r 2 (a,T 2 ) = J U'(p,a;k 1 ,k 2 ,k 3 ,T)w r 1 (p,T)diy(p;k 1 ,k 2 ,k 3 ,T), 

as an identity for generalized eigenvectors. Now taking inner products with k^ m on both sides 
gives (|4.2|) for the kernel of U' . So the kernel of our operator U is the same as the (formal) 
Racah coefficient. 

From here on we use the notation 

U{p,a;ki,k 2 ,k 3 ,T) = U^ a (T,p; ki,k 2 ,k 3 ). 

We are now ready to determine an explicit expression for the Racah coefficient. To write 
down the result it is convenient to introduce the following function: 

T{d + b + c + iX) 

(px[x;a,b,c,d) ~ 



r(o + b)T(a + c)r(l + a - d)T(l -d- iX)T(b + c + iX±ix) (4.4) 
x W(a + b + c — 1 + iX; a + ix,a — ix,a + i\,b + iX,c + iX), 

where Bailey's W notation is used for a very-well-poised 7-Fg-function, i.e. 

ttt / 1 j p\ ^ ( a, l + \a, b, c, d, e, / \ 

W(a;b,c,d,eJ) = 7 F 6 2 ;1 

\ 2 a ) l+o — 0, 1 + a — c, 1 + a — d, 1 + a — e, 1 + a — / / 

_ ^ a + 2n (a)n(b)n(c)n(d)n(e)n(f)n 

^ a n\(l + a - b) n (l + a - c) n (l + a - d) n (l + a - e) n (l + a - f) n 
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The parameters a, b, c, d are called the dual parameters, and they are defined by 

a = -(a + b + c + d- 1), b = -(a + b - c - d + 1), 

I . i ( 4 - 5 ) 

c= -(o-6 + c-d + l), d= -(a-b-c + d + 1). 

The function cp\(x;a,b,c,d) is called the Wilson function. Another useful expression for the 
Wilson function is obtained from transforming the yi^-function into a sum of two balanced 
4F3-functions, see ^1 (3.5)], 

a < t r\ T(l-a-d) 
4>x{x;a,b,c,d) - 



T{a + b)T(a + c)T(l -d± ix)T(l -d±i\) 
a + ix, a — ix, a + iX, a — i\ 



x 4^3 , , , , 

a + b,a + c,a + a 

r(a + d-i) ( 4 - 6 ) 



r(l + 6 - d)r(l + c - d)T(a ± ix)T(a ± iA) 

1 — <i + ix, 1 — (f — ix, 1 — <i + i A, 1 — (i — iX 
l + b-d, l + c-d,2-a-d 



X 4^3 



From this expression it follows that the Wilson function satisfies the duality property 

4>\(x;a,b,c,d) = (f> x (X;a,b,c,d). (4.7) 

Also we see from (j4.6|) that the Wilson function is invariant under a <-> 1 — d and b <-> c. For 
large x we have |12| 

d—a—b—c—2i\„irx(„ „,—2i\ , „ ^.2iA 



A (x; a, 6, c, d) ~ j^-^e" ^cix -AA + c 2 x 

where ci,C2 are independent of x. We show that the Racah coefficient U(p,a;ki,k2,k 3 ,r) is a 
multiple of a Wilson function. 

Theorem 4.5. T/ie Racah coefficient U(p,a;ki,k2,k 3 ,r) can explicitly be written as 

2n r(2k 1 )r(2k 2 )r{2k 3 )r(±2iT) , . , . , , 1 , 1 , 

) v ^ p \ </> p (a; h + it, k 2 + k 3 - -, k 3 - k 2 + -, 1 - h + tr). 

r(fe 2 - k\ - k 3 + ~ ± zr) ^ ^ 

For the proof we need the following lemma, which is Corollary 6.4 in |12| . 

Lemma 4.6. The Wilson function satisfies: 

4>\(x;a,b,c,d) = 

^ u f f-n,f + iX,f-iX \ ( -n, f - a, f + d - 1 \ 

u>/iere 

(/ + c-c±ix) n r(/±*A) 



C n (x, A) 



n! r(l + a - d + n)T(f + 6)r(/ + c)T(a ± iA)r(l - d ± i A) 



Proof of Theorem \4-5\ We use 1)4.2(1 . where we write the functions K r n m in terms of continuous 
dual Hahn polynomials and we put m = r = 0. Using the orthogonality relations for the 
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continuous dual Hahn polynomials s n (p; k 2 — k\ + |, k\ + k 2 ~~ 2' — ^ 2 + I) we obtain 

r(&3 + ir ± ip)T(k 3 — it ± ip)T(k 2 — ki — k 3 + | ± ir) 



= U(p,a;k l ,k 2 ,k 3 ,T)- 



2tt r(2fe 3 )r(A:2 - fci + i ± ip)r(±2ir) 



/ ( fc 2 ~ fc 3 + r + 2 ± »g)n f , . , . , , , lv 

> V W , , , s «n(r; fei + ^cr, fci - la, k 2 - h - k 3 + -) 

n\{2k 2 ) n 2 

x s n (p; fcx - k 2 + -, fei + fc 2 - -, fc 2 - h + -) 
(fcg ~ ^3 + 1 + icr) n (k 2 - k 3 + \ - ia) n 

n=0 




n\{2k x ) n 

+ -I 4- in. />'■> - 4- 4 — /'n ^ 



-n, k 2 - ki + \ + ip,k 2 - ki + \ - ip 



2fe,l 

n, fc 2 - ^1 - ^3 + i + ir, A; 2 - fci - fc 3 + i - ?t \ 
&2 - &3 + 5 + «cr, *2 - &3 + 5 - «c J ' 

provided that this function 1/ is an element of Ti(p}[(p; k\, k 2 ))- Here we assume for simplicity 
that r G [0, 00). In case r is in the discrete part of the support of the corresponding measure, 
the proof runs along the same lines. Now we apply Lemma 14.61 with parameters 

(a, b, c, d, f, x, A) h-> (k x + it, k 2 + k 3 - i, fc 3 - k 2 + ^, 1 - k x + ir, fc 2 - fci + ^, a, p), 

then we find the expression given in the theorem for the Racah coefficient. Using Stirling's 
asymptotic formula for the T-function and the asymptotic behaviour of the Wilson function, we 
verify that indeed V € TL{p\(p; k\, k 2 )) . Now the theorem follows from the fact that a function 
in H(pi(p; ki, fc 2 )) is uniquely determined by the coefficients in its expansion in the continuous 
dual Hahn polynomials s n (p; k 2 — k\ + 5, k\ + k 2 — 5, k\ — k 2 + D 

Remark 4.7. We assumed for simplicity that k\,k 2 , k 3 are such that discrete terms do not appear 
in the decompositions of 7r^. ® tt^, i = 1,3. If there are discrete terms in the decomposition, 
the Racah coefficients can be determined in exactly the same way, only discrete terms should 
be added to the integrals. In case the discrete terms correspond to positive discrete series, the 
Clebsch-Gordan coefficients for the tensor product of two positive discrete series are also needed 
here. 

4.2. The Wilson function transform of type I. Theorem 14.51 gives a natural interpretation 
of the Wilson functions in the representation theory of the Lie algebra su(l, 1). In ^2] the Wil- 
son functions are studied as eigenfunctions of a certain difference operator. Spectral analysis of 
the difference operator on a specific Hilbert space leads to a unitary integral transform that has 
the Wilson as a kernel. In this section we show how the integral transform, called the Wilson 
function transform of type I in |12j . is related to the interpretation of the Wilson functions as 
Racah coefficients. 

Let the parameters a, b, c, d £ C satisfy a, b, c, 1— d > 0, or b, c > 0, a = 1— d and K(a) > 0. We 
define the Hilbert space M. = M(a,b,c,d) to be the Hilbert space consisting of even functions 
that have finite norm with respect to the inner product 



(f,9)M= f(x)g(x)dm(x;a,b,c,d), 
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where m is the measure given by 



oo 







N , , , ,s 1 f°° Ma±ix)T{b±ix)T(c±ix)T{l-d±ix) 

f(x) dm{x; a, b,c,d) = — / f{x) dx. 

2tt J T(±2ix) 



We define an operator T = J 7 (a, b, c, d) on Ad by 

(JT/)(A) = (f,4> x (-,a,b,c,d)) M - 
This operator is initially defined on the domain consisting of functions / € A4(a, b, c, d) such that 
the integral (/, 4>\}m converges. We call this integral operator T the Wilson function transform 
of type I, or just the Wilson transform I. The following theorem is Theorem 4.12 in |12| . 

Theorem 4.8. The operator T [a ^b^c^d) : A4(a,b,c,d) — > A4(a, b, c, d) extends to a unitary 
operator, and the inverse is given by J- (a, b, c, d) . 

We sketch a proof of this theorem using the current interpration of the Wilson functions as 
Racah coefficients. 

Let us define the integral transform Z = X(ki, k 2 , ^3, r) by 

1 : C(u(p;ki,k 2 ,h,T)) -> C(u(a; k 3 , k 2 , h,r)) 



( I f)( (J )= / f(p)U(p,a;k 1 ,k 2 ,h,T)du(p;k 1 ,k 2 ,h,T). 
Jo 

Then from comparing squared norms of the funtions K^ m it follows from (|4.2|) that I is a partial 
isometry. By the following lemma we find that I extends to a unitary operator. 

Lemma 4.9. The set of functions {K^ m (p;ki,k 2 ,k3,T) \ n,m G Z>o, r G Z} spans the space 
C(u(p;k 1 ,k 2 ,k 3 ,T)). 

Proof. Let / G C(y{p; ki,k 2 , k 3 ,r)) be a function such that 



f{p) K n,m( T >P'i h,k 2 ,k 3 )dv(p; ki,k 2 ,k3,T) = 0, 
for all n, m, r. We integrate over all r, then by Fubini's theorem we have 

f(p) K Lm( T iP'i k 1 ,k 2 ,k 3 )dm r (T,p;k 1 ,k 2 ,k 3 ) = 0, 



for all n,m,r. Since the functions K^ m , n,m £ Z> form a basis for 7i{m r (x,y; ki, k 2 , k 3 )) we 
find that / = almost everywhere. □ 

Writing out the weight functions for the measure v{p; ki,k 3 ,k3,r), we see that this measure is, 
up to a multiplicative constant, the same as the measure m(p; a, b, c, d) for the Wilson transform 
I, where the parameters a, b, c, d (the Wilson parameters) are given by 

(a,b,c,d) = (fci + ir, k 2 + k 3 - - , k 3 - k 2 + -, 1 - fci + it). 

So the above defined integral transform I is the Wilson transform I. 

Note that the variables p and a came from decomposing the tensor product ir + <g> ir~ into 
a direct integral over the principal unitary series. So the support of the measure m for the 
Wilson transform I corresponds to the decomposition of ir + ®ir~ into irreducible representations. 
Discrete terms can appear in the decomposition, corresponding to discrete series or exactly 
one complementary series, and in this case the measure m for the Wilson transform I has 
corresponding discrete mass points. These are exactly the cases that one of the parameters 
a, b, c, 1 — d is non-positive. 

To define the Racah coefficients we might as well have started with the inverse of (|4.2I) 

i>oo 

K n,m( T - ,P',h,k 2 ,k 3 ) = / K^ n (T,a;k 3 ,k 2 ,ki)U(p,o-,ki,k 2 ,k 3 ,T)du(a;k 3 ,k 2 ,h,T). 
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Now interchanging k\ «-> k 3 , p <-> a, n <-> m, and comparing the obtained expression with (14,2(1 
gives the following property for the Racah coefficients: 

U(p,a;ki,k 2 ,k3,r) = U(a,p; k 3 , fc 2 , r). 

From l|4.5() we see that interchanging fci and /C3 is the same as interchanging the Wilson param- 
eters with the dual Wilson parameters. So the duality property (|4.7|) for the Wilson function 
can be obtained as a consequence of the associativity of the tensor product ® 71^ < X )7r fc 3 - This 
duality property can be considered as the main reason why the transform F is self-dual, up to 
an involution on the parameters. 



5. Racah coefficients for 7r + <g> ir p tg> ir 

In the same way as in the previous section we calculate the Racah coefficients for the tensor 
product representation tt^ <8> iTp £ <8> ft^ . Also we relate the generalized orthogonality relation 
for the Racah coefficients to the Wilson function transform of type II. 

5.1. Racah coefficients. The subspace of £ 2 (Z>q) (g) ^ 2 (Z) ® ^ 2 (Z>o) spanned by elements 
Cni 55 e n2 <8) e„3 for which n\ + n 2 — re 3 is constant, is invariant under the action of the Casimir 
in the tensor product Tvj^ <8> vr^ e ® 7r^. For this reason we define, for r£Z, 

jLr _ m+r ~ _ ~ jr-n 

rv n,m — i/n w c m — c " w l m 

Here and In are defined by (|H.4j) and Now the subspace 



A r = span c {/c nm \n,me Z>o} 

is invariant under the action of £1. 

We define corresponding functions K r n m by 

K, m (r, <r, ^p^,^) = (-l) m - r G^ +r (a-k 1 ,p,E)G- r (r; k 3 ,a,-h-e), (5.1) 
and we denote by m r (r, a; ki,p, e, £3) the corresponding orthogonality measure; 

f(T,a)dm r (T,a;k 1 ,p,e,k 3 ) = j j f(T,a)d(p^ +r (a;k 1 ,p,e) x p~ r (r; fe 3 , cr, -fei - e)). 

We should be careful with the definition of the function K T n m and the corresponding measure 
dm r when a is in the discrete part of the support of dp 2 (a; k%,p, e), since in this case the second 
function G~ r in (|5.1|) must be replace by a function F^ 3 , see also Remark 15.21 We leave the 
details to the reader. Now the following proposition is proved in the same way as Proposition 

ED 

Proposition 5.1. For r £ Z t/ie operator J\ : A r — > TC(rn(r, cr; fci, p, e, fc 3 )) defined by 

Ji : A£ 1 ^ ( (r,cr) i-> i££ (r, a; /ci, p, e, fc 3 )), 



is unitary and intertwines 7r^ ® vr^ e ® 7r fca ((A (g) 1) o A(f2)) luii/i M r2+ i . 
T/ie operator J 2 : A r — > 7i(m(T, C; &3, p, — £, ^3)) defined by 

J2 ■ K,m ^ (V>0 >-» ^m'n( T >C;^3,P,-e,^l)), 

is unitary and intertwines 71^ ® 7r^ E ® 7rj^((l ® A) o A(f2)) luif/i M t2+ i . 
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Remark 5.2. Using Propositions 13. 41 and !3.7l we can decompose the tensor product 7r^ (g)7r^ e (8)7r fe 
into irreducible representations in two ways. This gives 

/ oo 



oo / oo \ 

y* (y '^ e +fci-fc S dr ©© 7r fe-fci- e +j J dfj 




vr^ e+A . 1 _ fc3+j (iT©07r+ +A , i _ fc3+i _ / | . 



3 • j 
n V £ / 




©S / < fcl+e - fc3 - J -^©0 



^fcs-fci-e+j-J 







We see that the values for r corresponding to positive discrete series representations only occur 
for discrete values of a and for continuous values of Q (i.e. Q E [0, oo)). More precisely, if 
a = i{\ — k\ — e + j), then r can take the value r = i{\ — k\ + k% — e + I + j), where I € Z>o is 
such that t £ iR<o. A similar statement holds for values of r corresponding to negative discrete 
series representations. 

So we should be careful when writing down the function K r n defined by (|5.1j) for a = 
i{\ — ki — e + j), since then the second function G~ r in (|5.1|) must be replaced by a function 
'(t; fei + £ + /c 3 ), coming from the tensor product 7r^ + . ■ (g> 7rj~ . 

We first concentrate on the case r E [0, 00). We define the measure v t (<j; k%, p, e, /c 3 , r) to 
be the measure m r (T,a;ki,p,£,k 3 ) restricted to a line on which r = constant. For r — > 00 
this measure has an infinite number of discrete mass points (cf. Remark 13. 5|) . and we denote 
the corresponding measure by u, i.e. we omit the superscript 'r\ As in subsection 14.11 we can 
define the Racah coefficients as the kernel in the integral operator U : C(y(a; k\, p, e, fc 3 , r)) — > 
£(v((;fo,P, -e,ki,r)) mapping K^ m (r, a; ki, p, £, k 3 ) to K~^ n (r,C',h,P, s,ki). So the Racah 
coefficient U(a, (; ki, p, £, k 3 , r) is defined by 



K m r ,n( T , C k 3, P, -e, h) = J U (a, (; k u p, £, k 3 , r)K r n >m (r, cr; kx,p, £, k 3 ) d^(cj; fei, p, e, fe 3 , r). 

(5-2) 

In Remark 15.21 it is explained that discrete values of r only occur when one of the variables 
a, C is purely continuous and the other is purely discrete. In 1)5.2(1 the meaure dv is absolutely 
continuous (supported on [0, 00)) for r = i{\ — k\ + k 3 — £ + j) £ iM<o, and the measure dv is 
purely discrete (supported on the points 07 for Z E Z such that cj^ = £(5 — &i — e + Z) £ £M<o) for 
t = i(i + fei - Zc 3 + e + j) £ £K< . 

We use (|5.2(l and Lemma 14.61 to determine an explicit expression for the Racah coefficient 
as a multiple of a Wilson function. It turns out that for the discrete values of r the Wilson 



WILSON FUNCTION TRANSFORMS RELATED TO RACAH COEFFICIENTS 23 

function that we need can be simplified to a multiple of a Wilson polynomial 40 . The Wilson 
polynomials is a polynomial in x 2 which is defined by 

/— n, n + a + b + c + d — 1, a + ix, a — ix 
W n (x;a,b,c,d) = (a + b) n (a + c) n (a + d) n A F 3 [ , , , ;1 



a + b, a + c, a + d 

By Whipple's transformation [2j Theorem 3.3.3] this polynomial is symmetric in the parameters 
a, b, c, d. The following properties of the Wilson function turn out to be useful. 

Lemma 5.3. 

(i) The Wilson function has the following symmetry property: 

4>\(x; a + iu, b + iy,b — iy,l — a + ioj) = 4> u {y; a + i\,b + ix,b — ix,l — a + iX). 

(ii) The Wilson function 4>\{x; a, b, c, d) is symmetric in a, b, c, 1 — d. 

(iii) For A n = i(a + n), n G Z>o, the Wilson function is a multiple of a Wilson polynomial; 

, ( u rl\ = (-l) n W n (x;a,b, c, d) 

cp Xn \x, a, o, c, a) r(a + 6 + n)r(a + c + n)r(b + c + n)r(1 _ d±ix) 

Proof. The first property can be read off directly from the definition (|4.4j) of the Wilson function 
after the substitution (a, b, c, d) \— > (a + iu;, b + iy,b — iy,l — a + iu;). The second property is 
proved in [12| Remark 4.5(iii)]. The third property follows directly from (|4.6J1 . where for A = A rt 
the second 4.F 3 -function vanishes because of the factor T(d ± iA) -1 . □ 

Theorem 5.4. T/ie Racah coefficient U(o~,£;k±, p,£,k 3 ,r) can be expressed as a Wilson poly- 
nomial or a Wilson function as follows: 

• for t G [0, do), 

Ki 4>^(a; k\ + ip, k 3 + ir, k 3 - ir, 1 — k\ + ip), 

• /or r = i{\ - ki + £; 3 — s + j) G iM <0 and (" = i(| — k 3 + e + I) G iM<o, 

^2 Wj^i(a; 2k 3 - k\ - - - e - j, h + ip, k x - ip, - - h - e - j), 

• for t = i(^ + k\ — k 3 + e + j) G iM<o and <r = i(^ — k\ — e + I) G iM<o, 

#3 Wi-KC; 2fci - fc 3 - - + e - j, k 3 + ip, k 3 - ip, ^ - k 3 + e - j), 

u>/iere 

i^i =2sinvr(i - iti - e + k 3 ± iT)T(2k 1 )T(2k 3 )T(±2iT)T(^ - e + k 3 ±iQT{]^ + k x + e ±ia), 

_ r(2fc 3 - 2fci - 2e - \)Y(2k 3 )Y(2k 3 -2s -I) (— l) J "~ f Z!(^ - fc 3 + fci + e) 3 - 

2 r(2A;i - 2A: 3 + 2e)r(2A: 3 -\-e-l±ip) (§ - fc 3 + A* + e) 3 (± + A* + e ± ia)^{2k{)~i ' 

_ T(2fci - 2fc 3 + 2e - l)r(2fci)r(2fci + 2e-l) (— l) J "~ f Z!(^ - fci + fc 3 - e) 3 

M T(2A; 3 - 2fci - 2e)r(2fci - | + e - I ± ip) (f - h + k 3 - e) 3 -(i + A; 3 - e ± icr)j(2fe 3 )j_ z 

Proof. First assume that r G [0, oo). From ()5.2|) with m = r = and using the orthogonality of 
the continuous dual Hahn polynomial s n (a; k\ + ip, k\ — ip, \ — £ — k\), we find 

Uc_{a;ki,p,e,k 3 ,T) = 

r ~ £ ± jp)n F / -n, 5 -e-fci +io-, | -e-A;i -icr \ 

^ n!(2fe)„ 3 2 ^ i- £ + ip,i- £ -ip ; J 



3^2 



-n, I — e — fei + A; 3 + ir, | — e — k\ + k 3 — ir 
\ - e + k 3 + i(, 4 - e + fe 3 - i( 
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Using Lemma 14.61 with parameters 

(a, b, c, d, f, x, A) i-> (ki + icr, k 3 + i£, k3 — — k± + ia,^ — e — k\ + k 3 ,p, r) 

we find that the Racah coefficient is a multiple of the Wilson function <f> T (p; k\ + ia, k^ + ?C> k% — 
i£, 1 — ki+io). Then the theorem follows from applying Lemma [5.3f i). and writing out explicitly 
the factor C. 

For discrete values of r the proof runs along the same lines, using also Lemma l5.3f ii) and 
(iii). □ 

Remark 5.5. Note that the Wilson function and the Wilson polynomials in the theorem do not 
depend on e. 

5.2. The Wilson function transform of type II. The Wilson function transform of type 
II |12j is an integral transform depending on the four Wilson parameters a, b, c, d and an extra 
parameter t. Let the parameters a, b, c, d, t satisfy 

9fc(a),9fc(6) > 0, a = l-d, b = c, t £ R. 

For these parameters we define the Hilbert space H = 7i(a, b, c, d; t) to be the space consisting 
of even functions that have finite norm with respect to the inner product 



(/, 9)u = J f(x)g(x) dh(x; a, b, c, d; t), 
where h is the measure defined by 

/s~t poo 
fix) dh(x; a, b, c,d;t) = — / f(x) H(x) dx + iCS" f(x) Res H{z 
27r Jo z=x 



xet> 



where 

r(a ± ix)T(b ± ix)T{c ± ix)T(l — d ± ix) 



H(x) = H(x; a, b, c, d; t) 



smir(t ± ix)T(±2ix) 



C = C(a, b, c, d; t) = \J sin-7r(a + t) sin7r(6 + t) sin7r(c + t) sin7r(l — d + t), 
and T> is the infinite discrete set 

V = {i(t - n) | n G Z, t - n < 0}. 
We define the dual parameter t by 

t = l-b-c-t, 

and the dual parameters a,b,c,d are still defined by (|4.5|) . Now we define an operator Q = 
Q(a, b, c, d; t) by 

(Gf)(\) = {f,M-,a,b,c,d)) H , 

for all functions / S Ti such that the integral on the right hand side converges. The operator Q 
is called the Wilson transform of type II, or just the Wilson transform II. The following theorem 
is Theorem 5.10 in ^5] 

Theorem 5.6. The operator Q : Tl(a, b, c, d; t) — > 7i{a, b, c, d; t) extends to a unitary operator, 
and its inverse is given by Q{a,b,c,d;t). 

We will sketch how the unitarity of Q can be obtained from the interpretation of the Wilson 
functions as Racah coefficients for the tensor product 7^ <8> 7r^ £ <S> vr^ . 

Let r € [0, 00). Let us define an operator 2 = I(k\, p,e,ks,r) : C(v(a; k\, p, e, ks, r)) — > 
^{^(C)h,P,-£,ki,r)) by 



( I f)(0= / f(v)U{o-,(;k 1 ,p,e,k 3 )dis(a;k 1 ,p,E,k 3 ,T). 
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Then, in the same way as in subsection 14.21 we can obtain from (|5.2|) that X is unitary, and 
its inverse is the same tranform after k\ <-> k% and e » —e. Comparing the measure v with 
the measure h for the Wilson transform II, and using Theorem 15.41 we see that we have, for an 
appropriate function /, 



where 



and 



(X/XC) = / </>d a '> a ' & ' c ' d )A(a)f(a) dh(a; a, b, c, d; t), 



(a, 6, c, d, t) = (ki + ip, k 3 + ir, /c 3 - ir,l — kx + ip, - 



A{a) =T(~ + fci + e ± iff)" 1 , 5(C) = r(i - e + fc 3 ± iQ~\ 



E 



T(\-k z + e + k l ± %t)V{\ +e±ip) 



T(±-k 1 -e + k 3 ± iT)T(\ -e±ip) 

So we have = £7 {^b(0 ° & ° ^4(-)/)(C)i where <7 denotes the Wilson transform II with 

the above parameters. Note that A(-) = B(-) and E = E . Here we use the notation O = 
0(d, b, c, d, i) for any object O = 0(a, b, c, d, t) depending on the parameters a, b, c, d, t. Observe 
here that replacing (a,b,c,d,t) by their dual parameters is the same as replacing (ki,k 3 ,e) by 
(&3,fci,— £■). The operator M^(.) is invertible, because -A(c) ^ for a in the support of the 
measure v{o\ k±,p, e, &3, t). 

The measure ^ = u(-; k\,p, e, k 3 , r) is related to the measure h(-; a, b, c, d; t) by 

f{a)dv(a)=N j f(a)\A(a)\ 2 dh(a;a,b,c,d;t), 

where, if C is the normalizing constant in the definition of the measure h satisfying C = C, 

T(± + ki + e - k 3 ± ir) 
~ 2CT(±2ir)r(2A; 1 )r(2A;3)r(i -e±ip)' 

So we have f2)c{u) = iV(-4 /1, -4 f2)H- On the other hand, we find from the unitarity of I, 

{h,h)c.{y) = @hiZh)c{y) 

J! \ A-l r>f A J- \\ 



| J e| 2 (I- 1 g(^/ 1 ),l- 1 g(A/ 2 )) £(p) 

|£| 2 JV(0(A/i),S(A/ 2 )>, 



Now from |_E| 2 iV = JV it follows that the Wilson function transform Q is unitary. 

Let us also remark that the infinite set of discrete mass points of the measure h for the Wilson 
transform II corresponds to the infinite sum of discrete series in the tensor product decomposi- 
tion of -k^ ®tt p . The absolutely continuous part of h corresponds to the direct integral over the 
principal unitary series. 

For r [0, 00) the Wilson functions in Theorem l5.4l reduce to Wilson polynomials. The Wilson 
polynomials are orthogonal with respect to a positive measure jlU]. Let the parameters a, b, c, d 
be such that that non-real parameters appear in pairs of complex conjugates with positive real 
part, and such that the pairwise sum of the real parameters is positive. Then the orthogonality 
relations for the Wilson polynomials read 



00 



j W n (x)W m (x)w(x)dx + W n (xj)W m (xj)wj = h n 6 nm . 



.1 
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where we denote W n (x) = W n (x; a, b, c, d). The sum is over all j € Z>o such that Xj = i(e + j) € 
iM<o, where e is the smallest of the real parameters. In particular, the sum is empty when there 
is no negative parameter. The weights and the squared norm are given by 

r(a ± ix)T(b ± ix)T(c ± ix)T{d ± ix) 



w[x 



Wj =Res w(x 



x=x 



T(±2ix) 

_T(a + b + n)T(a + c + n)T(a + d + n)T(b + c + n)T(b + d + n)T(c + d + n) 
n ~ r(a + b + c + d + 2n) 

x (n + a + 6 + c + <i — l) n n!. 

The Wilson polynomials form an orthogonal basis for the Hilbert space consisting of even func- 
tions that have finite norm with respect to the inner product 

(f,g) = 7T f{x)g{x)w(x)dx + y^f{xj)g{xj)wj. 

So we have the following pair of unitary integral transforms 

\Wf)(n) = (f,W n ), 

oo 

f(x) = Y / (Wf)(n)W n (x)h- 1 

n=0 

We call W = W(a, b, c, d) the polynomial Wilson transform. 

It can be checked that the unitarity of polynomial Wilson transform in the case a, d > and 
b = c can be obtained from the interpretation of the Wilson polynomials as Racah coefficients in 
Theorem l5.4l We leave this to the interested reader. Another interesting fact that we obtain from 
this interpretation is that the polynomial Wilson transform W(a, b, c, 1 — d) can be considered 
as a degenerate case of the Wilson function transform of type II G(a, b, c, d; t), namely the case 
that the parameters satisfy 

a = l-d, K(a) > 0, b,c£R, b + c>0, t = n-c+l, 

for n € Z>o- 

6. More Racah coefficients 

In this section we determine some more Racah coefficients for 3-fold tensor product represen- 
tations of su(l, 1). 

Let us start with the tensor product representation <£> ir^ <g> ir^ . We assume for simplicity 
that there are no discrete terms in the decomposition of 7r^ ® ir^ . The ^-invariant subspace in 
this case is 



A r = span c {e ni ® e n2 <g> e n3 | m, n 2 , n 3 G Z> , n\ + n 2 - n 3 = r} C l 2 (Z>o)® 3 . 

Using the results from section |21 we obtain operators J\ and J2 as in Propositions 14.11 and 14.21 

Let us define the measure mi by 

fk fk OO n 

If f(T,j)dm 1 (T,j;k 1 ,k 2 ,k 3 ) = Y] / f(r,j)W{j;ki,k 2 )dfii(T;k 1 +k2+j,h), 

and let dm^ be the product measure defined by 

// f(r,p)dm 2 (T,p;k l7 k2,k 3 )= // /(r, p) d[m(z; k 2 , k 3 ) x fj, 2 (x; k 1: p, k 2 - k 3 )). 
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Now the operators J\ and J 2 are given by 

Ji : A r -> n(m 1 (T,j;k 1 ,k 2 ,k 3 )) 
e ni <8>e„ 2 <g> e„ 3 ^ Mr, j) i-> <J +r (i; fa, fc 2 )F^ J (r; fcj + fa+ j,k 3 )^, 



and 



J 2 : ^ r -> U{m 2 (r,p; fa, k 2 , k 3 )) 
2ni ®e„ 2 (8>e n3 1 ^ ((r,p) i-> F^~ ni (p;k 2 ,k 3 )G r ni (T;fa, p,k 2 - k 3 ) 



These operators intertwine yrj^ ®7r^ ® 7r fc 3 ((^ ® 1) A(f2)), respectively ® 7T^ <g>7rjjT ((1 <8> A) o 
A(f2)), with M r 2 + i . We restrict the measure mi(x, y; fa, k 2 , k 3 ) to a line on which r = constant, 
then we obtain a measure which is supported on Z>o- We will denote a mass point of this 
measure at j € Z>o by w(j; fa, k 2 , k 3 ,r). Furthermore, we define the measure v(p; fa, k 2 , k 3 ,r) 
to be the measure m 2 {r, p; fa,k 2 , k 3 ) restricted to a line on which r = constant. Now the Racah 
coefficients U(j, p;fa,k 2 ,k 3 ,r) can be determined from 



(p; k 2 , k 3 )G r ni (r; fa, p, k 2 - k 3 ) = 

^2u(j, p; fa,k 2 , k 3 ,T)v^ +r (j; fa, k 2 )F,Jl~ J (T; fa + k 2 + j, k 3 ) w{j; fa,k 2 , k 3 ,r), 6-1 

or equivalently from 

<? +r (i; fa,k 2 )F^{r; fa + k 2 + j, k 3 ) = 







U(j,p; fa,k 2 ,k 3 ,r)F^ 3 ni (p; k 2 ,k 3 )G r ni (T;fa, p,k 2 - k 3 )dv(p; fa,k 2 ,k 3 , 



(6.2) 



T . 



Then from Lemmas 14.61 and 15.31 we find the following expression for the Racah coefficient. 
Theorem 6.1. The Racah coefficient U(j, p;fa,k 2 ,k 3 ,r) is the Wilson polynomial 

K Wj(p; k 2 -k 3 + ^,k 2 + k 3 -^,fa + ir,fa- ir), 

where 

2 (-iyT(±2iT)T(2fa + 2k 2 + 2j)T(2k 3 ) simr(fa + k 2 - k 3 + \ ± ir) 



K 



(2fa)jT(fa + k 2 + k 3 - \ + j ±ir)J (fa + k 2 - k 3 + \ ± ir)j(2ki + 2k 2 + j)j 



Proof. We use formula (|6.2|) , We use the orthogonality relations for F^~ ni (p) to expand the 
Racah coefficient in terms of the Clebsch-Gordan coefficients. Then by putting m = and 
r = j and writing out the Clebsch-Gordan coefficients explicitly, we obtain 

U(j,p;fa,k 2 ,k 3 ,T) = 

c , (j + l)n 3 (2fc 2 + j)n 3 F f -«3, fa - k 3 + j + \ + ip,k 2 - k 3 + j + \ - ip 
^ n 3 \(2k 3 ) m 3 2 V 2k 2 +j,j + l 

"3 — U 

/ -n 3 , fa + k 2 - k 3 + j + \ + ir,fa + k 2 - k 3 + j + \ - ir \ 
X 3 r 2 I ; J- . 

V l,2fa+2k 2 + 2j J 

Applying Lemma 14.61 and Lemma l5.3l i) shows that the Racah coefficient is equal to the Wilson 
function C <p i ^ kl+k , 2 _ij r ^{p', fa + k 3 — ^, fa + ir, fa — ir, k 3 — k 2 + ^), where the factor C is given 
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by 

2r(±2ir)r(2A; 1 )r(2A:3)r(A:3 - k 2 + - ± ip) 

T(2k 1 + 2k 2 + 2j) sin ir(ki + k 2 -k 3 + \±iT) 
x ■ 



'(fci + k 2 - k 3 + \ ± iT)j(2ki + 2& 2 + j^- 
Finally the theorem follows from applying Lemma 15.3( 11) . □ 



Remark 6.2. Koornwinder [22] found a similar group theoretic interpretation of the Wilson 
polynomials by establishing the Wilson polynomials as connection coefficients between two 
0(p) x 0(q) x 0(r)-invariant bases on a hyperboloid in M p+!?+r . 

Finally we give the Racah coefficients for the tensor product representations 7r^ <8> 7rf e (£> tt^ 
and 7r^~ <8>7r^ 2 ® 7r ^ e - The derivation of these Racah coefficients is completely the same as before, 
therefore we omit all details. 

The Racah coefficient U (a, £; e, £3) for the tensor product 7r^ (g) vr^ e (g) -zr^ is the kernel 
in the transform 2 : C(y{a; k%,p, e, kz, r)j — > £(z/(£; k 3 , p, e, &i,t)) , which is defined by 

(2 [a -» (G r n - n * (a; h ,p, e) G r n3 (r ; k 3 , a, e + h )] ) «) = 

G;;" 1 (C;A ; 3,P I £)G ! ; i (r;fc 1 ,C,e + A ; 3). 

Here the measure ^(<r; fei, p, e, &3, r) is the restriction of the product measure 

P2(cr;k 1 ,p,e,k 3 ) x p 2 (r;k 3 ,a,e + h) 

to a line on which r = constant. The Racah coefficients can be expressed as a multiple of a 
Wilson function 



U(a,(;h,p,£,k 3 ) 

rC2ki )r(2ko)r(l - F + i n \ 

-<^(<r; k 3 + it, h + ip, fci - ip,l - k 3 + it) 



2TrT(±2iT)T(2ki)T{2k 3 )T{\ -e±ip) 



V(\-k x -k 3 -£± it) 

The T-functions in the measure v depending on e are cancelled by the same T-functions in the 
expression of the Racah coefficient, so the corresponding integral transform 2 does not depend 
on e. This integral transform 2 is equivalent to the Wilson transform I. 

The Racah coefficient U(j, a; k\, k 2 , p, e) for the tensor product ir^ <g> 71^ ® tt^ £ is the kernel 
in the transform 2 : ^ 2 (Z>o, w(j; k\ , k 2 , p, e)) — > C(y(a; k\, k 2 , p, e, r)), which is defined by 

-> v^fthMW-is-jfcki + k 2 +j, p,e)])(o) = 

G^Znl-nsi^ k 2,P, s)G r ni (r; k 1} a, k 2 + e). 

Here w(j\ k\, k 2 ,p, e, r) is the mass at the points j £ Z>o of the measure 

W(j;k 1: k 2 ) x p 2 (T;h + k 2 + j,p,e) 

restricted to a line on which r = constant, and z^(cr; fcj, &2, p, e, t) is the restriction of the product 
measure 

p 2 (a;k 2 ,p,e) x p 2 (r,ki,a,k 2 +e). 
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The Racah coefficient can be expressed as a multiple of a Wilson polynomial, 

U(j,a;ki,k 2 ,p,e) = 

(-iy 2vr r(±2tr)r(| - £ ± ip)r(2Jfei + 2k 2 + 2j) 



(2fc 1 ) 3 T(i -k 1 -k 2 -e± tr)r(fci + fc 2 + j + tp ± ir)r(fci + fc 2 + j - ip ± it) 

— I — /=" — I — -i — I — 7.t\ „■ 

Wj(<t; /c 2 + ip, fc 2 - ip, fci + it, k\ - it). 



{k 1 + k 2 + e + \± ir) 3 
{2k 1 + 2k 2 +j) j " ,1 



The corresponding integral transform 2 does not depend on e, and X is equivalent to the poly- 
nomial Wilson transform. 



7. Integral and summation formulas 

7.1. Formulas related to 3- fold tensor products. In the previous sections we found ex- 
plicit formulas for the Racah coefficients by expanding the Racah coefficients in terms of the 
Clebsch-Gordan coefficients. Writing the Clebsch-Gordan coefficients as (continuous) dual Hahn 
polynomials, we obtain several summation formulas for these polynomials. 

Let us start with the tensor product ir^ <g> tt^ ® n^ 3 from section El We write the Clebsch- 
Gordan coefficients as (continuous) dual Hahn polynomials and the Racah coefficient as a Wilson 
polynomial, then formula (|6.1|) gives a summation formula involving Wilson, dual Hahn and 
continuous dual Hahn polynomials. Note that the sum is finite, since the dual Hahn polynomial 
T n (j; 7, 5, N) equals zero for j > N. Using the orthogonality of the dual Hahn polynomials, we 
obtain a summation formula which can be considered as a generalized convolution identity for 
continuous dual Hahn polynomials. 

Theorem 7.1. The continuous dual Hahn polynomials satisfy the following generalized convo- 
lution identity: 

n+j 



3 )(2k l ) j T^ki-l^hz-hn + j) 

l—n V / 



x S n+j -i(p; k 2 - k 3 + -,k 2 + h - -,k 3 - k 2 + I + r + -) 
x Si(t; k\ + ip, k\ - ip,k 3 - k 2 - k\ + r + -) 



1 , , 1 

-,k 2 + k 3 --, 

1 . , , I 



Wj(p; k 2 - k 3 + -, k 2 + k 3 - -, k\ + ir, k\ - it 



x S n (T; ki + k 2 - k 3 + j + -,fa + k 2 + k 3 + j - -, k 3 - k 2 - k\ + r + -). 

Remark 7.2. Several formulas with a similar structure, involving other hypergeometric orthogo- 
nal polynomials, are found by Koelink and Van der Jeugt |25j . These formulas and Theorem 17. II 
can be obtained as limit cases of one "master formula", a generalized convolution identity for 
Wilson polynomials, which is obtained by Lievens and Van der Jeugt |3L)l Thm.l]. In Theorem 
17.61 we show that the "master formula" itself also has an interpretation in the representation 
theory of su(l, 1). 

Next we consider the tensor product 7r^ <g> ir^ <S> 7r^ from section I1J Writing the continuous 
dual Hahn polynomials in (|4.2|) in their ususal normalisation ()3.2j) . and using Theorem 14.51 
and orthogonality relations for the continuous dual Hahn polynomials, we obtain an expansion 
formula for the Wilson funtions. 
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Theorem 7.3. The continuous dual Hahn polynomials and the Wilson functions satisfy 
<Mp; h + ir, h + k 2 - -, h - k 2 + -, 1 - k 3 + ir) 

xS m (r; k 3 + ip, k 3 - ip, k 2 - h - k 3 + r + -) = 

oo 1 i 1 

}Cn S n (p; ki - k 2 + - , h + k 2 - -, k 2 - ki + m + r + -) 

n=0 

x £ n (r; ki + id, /ci - ia, k 2 - ki - k 3 + r + ^) 

11 1 

x 5 m (a; fc 3 - k 2 + -, fc 2 + h - 77, fc 2 - k 3 + n + r + -), 



where 



2 7 2' 2 

r(fc 2 - ki + m + r + \ ± ip) 



n! r(2&i + n)r(2/c2 + n + r + n)T(n + r + m + l)r(/c3 — ir ± ip)T(k 3 + ir ± ip) 



The tensor product (g> 7T^ e (g> 7r fca leads to a summation formula which is equivalent to the 
formula in Theorem 17.31 The tensor product 7r^ ® tt^ 2 <g> 7r^ leads to a generalized convolution 
identity for (dual) Hahn polynomials, see |25| . 

Since the coproduct A is cocommutative we have tt\ (g> ir 2 = ir 2 ni for representations 
■Ki, i = 1,2. We use this symmetry of the tensor product to find relations between Racah 
coefficients. In order to make the symmetries in the tensor products more transparant, we will 
denote tti <g> tt 2 <g> tt 3 o (A <g> 1) o A by [wi ® tt 2 ] (g) tt 3 , and similarly for the operator related to 
(1 <S> A) o A. In order to distinguish between Racah coefficients for different tensor products we 
add superscripts +, — or P to the representation labels. For instance, U (a, (; kf, (p, e) p , k 3 , r) 
is the Racah coefficient for the tensor product representation n£ <g> ir P E ir^. We will use 
similar notation for the corresponding measures. Morover, we denote the corresponding integral 
transforms by Z([k± , (p, e) p ], fcT ,r) and Z(kf, [(p, e) p , k^\,r), which are each others inverse. 

From symmetries in the tensor product 7T^ <S> ir^ 2 ® ir^ we find the following theorem. 

Theorem 7.4. For k\ + k 3 — | > and k\ — k 3 + | > we have 
I roc J 1 

— / Wj(a; ki - k 3 + -, ki + k 3 - -, k 2 + ir, & 2 - »r) 
2vr Jo 2 2 



r(A: 2 + ir ± icr)r(A; 2 — ir ± ier)r(fei + k 3 — \ ± icj)r(/ci — k 3 + ^ ± id) 



x </> P (c; &2 + ir, h + k 3 - -, k x - k 3 + -, 1 - k 2 + ir 



r(±2ia) 

(— 1) j 'Wj(/j; k 2 - k 3 + /c 2 + /c 3 - ^, fci + ir, % - ir). 



f/(7 



r " " 2 

Proof. We can factorize the transform I(kf, [k^~, k 3 ],r) using the following chain of equivalent 
tensor products: 



7T 



I 



fc 2 ^ L"fci ^ "fc 3 J 



The T above an arrow indicates that the arrow corresponds to an integral transform with a 
Racah coefficient as a kernel. This gives us the following factorization, 

l(k+, [k+,k^],r) = Cl(k+, [k+, fcg ],r) ol(k+, r), (7.1) 
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for a (still to be determined) non-zero factor C. This factor can be found from applying both 
sides to the same function. We denote J = [kf , k 3 ], r) o X{k^ , [A;^, k% ], r), then from 

applying (j4.2j) and (j(12j) we obtain 

(J[P -» ^"HP); fe, fojG^ (r; A^p, fc 2 - fc 3 )]) (j) = v r + n *{j; k 2 , h)F^(r; h + k 2 + j, k 3 ), 

where r = n\ + ri2 — 723. Also, from (|6,2j) we obtain for T = Z(kf, [k£ , k 3 ], r) 

I[p - F„V ni (p); fca, fc 3 )G; a (r; p, k 2 - k 3 )] ) (j) = k u k 2 )F^(r; h + k 2 + j, k 3 ). 



ll +n Hr,h,k 2 ) = (-ly^iv^^k^h), 



Now C is obtained from 



which follows directly from (|2.6|) . 

We apply both sides of (|7.1j) to a suitable function /. If we choose / smooth enough we may 
interchange the order of integration on the right hand side, and since / was choosen arbitrary 
we then obtain 



(-iy^-U(j,p;kf,k+,k^r) 



U(j,a;k£,kf,k 3 ,r)U(p,a;kf ,k 3 ,k£ ,r) dv(a;k£ ,k 3 ,kf,r), 



where we use 



/ g{o)dv(a;k^ ,k 3 ,k£,r) = g(o) dv(a;kf ,k£ ,k 3 ,r). 
Jo Jo 

Writing the Racah coefficients as Wilson functions or Wilson polynomials, we obtain the theo- 
rem. □ 

We remark that Theorem 17.41 is valid for all ki,k 3 > if appropriate discrete terms are 
added to the integral, corresponding to the discrete terms in the tensor product decomposition 

Let us relabel 

(a,b,c,d) = (k 2 + ir, h + k 3 - -,ki - k 3 - -, 1 -k 2 + ir), 

then Theorem 17.41 can be written as 

(FWji- ; a, b, c, 1 - d)) (p) = (-iyWj(p; a, b,c,l- d). 

This is j!21 Theorem 6.7], which states that the Wilson transform I maps a basis of orthogonal 
polynomials in Ai(a, b, c, d) to a basis of orthogonal polynomials in M(d, b, c, d). 



7.2. Formulas related to 4-fold tensor products. In general the Racah coefficients, or 6j- 
symbols, are the recoupling coefficients that connect the coupled basis vectors for {tt\ ®tt 2 ) ®tt 3 
with the coupled basis vectors for n\ ® (tt 2 ®tt 3 ). The 9j-symbols are the recoupling coefficients 
related to a 4-fold tensor product. These 9j-symbols can be expressed in terms of the 6j- 
symbols. In this section we consider 9j'-symbols that connect the coupled basis vectors for 
[(tti <S> n 2 ) <S> 7T3] VT4 with the coupled basis vectors for tti (g> [ir 2 <S> (n 3 <8> ^4)]. There are several 
possible transitions from the first coupled basis to the second coupled bases. We consider the 
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following two transitions: 

X X 
[(iTi (g) 7T 2 ) gl 7T 3 ] (8) 7T 4 — > (7Ti g> 7T 2 ) (g [7T 3 ® 7T 4 ] — * 7Tl (g (7T 2 <g [7T 3 <g 7T 4 ]), 

[(7Tl (g> 7T 2 ) <g> VT 3 ] (g 7T 4 [7T1 <g> (7T 2 (8) 7T 3 )] (g 7T 4 (7.2) 

X X 

— * 7T1 <8) [(7T2 ® 7T 3 ) (8) 7T 4 ] — » 7Ti (g [7T 2 g> (7T 3 (g 7T 4 )]. 

These possibilities give rise to two different expressions for the 9j-symbols, so we obtain a 
summation or integral formula relating different Racah coefficients. Biedenharn j§] and Elliott 
|1U| used a similar method to find the famous Biedenharn-Elliott identity, or pentagon identity, 
for the su(2)-Racah coefficients. 

Let us consider the tensor product 7r^ (g> 71"^ (g 71"^ g> -k^ . The two transitions from (|7.2I) lead 
to the following identity for operators: 

l([kf,k+],(a,k 4 -k 3 ) p ,T) ol([(k 1 + k 2 +j) + ,k^},k+,r) 

= Cl([k£, k 3 ], fc+ C) o 2*([fc+ ( Wj A; 2 - fc 3 ) p ], fc+ , r) o I([fc+ fc+], fc 3 -, p) , 

where the factor C is still to be determined. 

This factor can be determined in the same way as in the proof of Theorem 17.41 Then we find 
that the factor we are looking for can be determined from 

CGT~\t; fc 4 ,p, ki+k 2 +j - k 3 ) = G r n (r; fc 4 ,p, k x + k 2 - k 3 ). 

and this gives 



C 



(ki + k 2 -k 3 + \± ip)j 



(ki + k 2 - k 3 + k A + \ ± ir)j 
Then, similar as in the proof of Theorem I7.4| we obtain the following identity. 
Theorem 7.5. For (, p,a,T G R, k\,k 2 ,k 3j k 4 > 0, j G Z>o, & 2 + /c 3 — | > 0, and /c 2 — fc 3 + | > 0, 



f°° 11 
/ 5(cj)Wj(a;; k 2 - k 3 + -,k 2 + k 3 - -,k 1 + ip, k x 
Jo ll 



ip) 



x (f) a (uj; fc 4 + i(, k 2 + k 3 - -,k 2 - k 3 + -, 1 - k A + iC) 
x 4> t {uj; k\ + ip, fc 4 + i£, /c 4 — i£, 1 — &i + ip)du 
Wj(C; k 2 + id, k 2 — ia, k\ + ir, k\ — it) 

x 4> a (p; k 4 + ir, fci + k 2 + fc 3 + j - -, fci + A; 2 - fc 3 + j + -, 1 - k 4 + ir), 



where 



. . _ i r(fc 2 + fe 3 - | ± icj)r(/c 2 - fc 3 + ^ ± zw)r(^i + ip ± iw) 

~ 2vr r(fci + £: 2 + /c 3 + j - ± ± zp)r(fci + fe 2 - & 3 + J + 5 ± ip) 
r(fei — ip ± iw)r(A: 4 + iC ± iw)r(&; 4 — i£ ± ioJ) 
X r(±2iw) ' 

Note that we used the symmetry relation Lemma l5.3lT j for the Wilson functions. If discrete 
mass points with the appropriate residues at the poles of the weight B{u) are added to the 
integral, Theorem 17.51 is valid for k 2 ,k 3 > 0. 
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Using the orthogonality relations for the Wilson polynomials, it follows that Theorem 17.51 is 
equivalent to the following bilinear summation formula for the Wilson polynomials: 

y^2fa + 2k 2 + 2j - 1 (2fa + 2k 2 - l)j 



where 



_ n 2fa + 2fa-l j\{2k 1 ) j (2k 2 ) j 

x (f) a (p; h A + ir, ki + k 2 + fa + j - -, fa + k 2 - fa + j + -, 1 - fa + ir) 

x Wj(uj; fa-fa + ^,k 2 + fa-^,fa + ip, fa - ip) 
x Wj(C; k 2 + ia, k 2 — ia, fa + ir, fa — ir) 

=K<f> a (uj; fa + i(, fa + fa-^,k 2 -fa + ^,l-fa + iQ 
x 4> t (lu; fa + ip, k A + i(, fa — i(, 1 — fa + ip), 

Y{2fa)Y{2k 2 )T{fa + i(± ico)Y(fa - i( ± iu) 



K 



Y{2ki + 2k 2 ) 

This sum can be considered as the q = 1 analogue of the bilinear summation formula for Askey- 
Wilson polynomials |26| Theorem 3.3] which is obtained by a different method. 

In the same way as we obtained Theorem 17.51 we can find Biedenharn-Elliott identities for 
other tensor products. From 7r^ (8 7rj^ (8 7r^ (8 ttT^ we obtain a generalized convolution identity 
for the Wilson polynomials, also involving Racah polynomials. 

Theorem 7.6. Let fa, fa, fa, fa > 0, ji,j 2 £ Z>o and p,Q £ 1, then the Wilson polynomials 
satisfy the following generalized convolution identity 

jl+j'2 

c j3 R j3 till 2k 2 - 1, 2fa - 1, ~h ~ 32 - 1, 2fci + 2/c 2 + ji + 32 ~ 1) 
^,fa + fe 4 - 

1 . . , 1 



x Wj 3 (/o; fc 3 - & 4 + -, & 3 + fc 4 - -, & 2 + i£, &2 - «C) 



where 



x Wji+j2-j 3 (C; ^2 + &3 - fe 4 + i3 + 2 > fc 2 + £3 + &4 + J3 - h + ir, fei - ir) 

: Wj 2 (p; A; 3 - fa + -, fc 3 + fa - -, fci + fc 2 + ji + ir, fa + k 2 + j x - ir) 
x Wj 1 ((;fa + ip, fa - ip, fa + ir, fa - ir) , 

31 + h\ (2fc 2 )j 1 (2fc 2 + 2fa) 2j3 (2fa + 2fa + 2fa + ji + j 2 - % 
J3 / (2fc3) i 3(i3 + 2fc2 + 2fc3-l)i 3 (2fc2 + 2fc3) il+i2+i 3 



Remark 7.7. This convolution identity is obtained in j3U] by analytic continuation of the classical 
Biedenharn-Elliott identity for su(2) Racah coefficients, which can be considered as a generalized 
convolution identity for the Racah polynomials. Many generalized convolution identities for 
hypergeometric orthogonal polynomials can be obtained from Theorem 17.61 as limit cases. For 
instance, if we put ir = fa + fa + fa — fa — ^ in Theorem 17.61 and we let fa — ► oo, then we obtain 
a generalized convolution identity for continuous dual Hahn polynomials which is equivalent to 
Theorem 17. II See [201 an d ES] for other limit cases. 

The Biedenharn-Elliott identity for the tensor product ir^ (8 ir~£ 2 (8 7r^ 3 (8) r^ leads to the 
following theorem. 
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Theorem 7.8. For p,(,T G E, &i, fe, fc 3 , fe* > 0, ji,j 2 € ^>o> an d k 2 + k 3 — \ > 0, 

/•oo 11 

/ B(u) W h (uj; h - k 2 + -,k 2 + h - -,k 4 + i(,ki - iQ 
Jo 22 

x Wj^uj; k 2 -k 3 + -, fc 2 + ^3 - 2' fc i + fc i ~ i P) 
x t (cj; fci + ip, ki + i(, k^ — i(, 1 — k% + ip) du; 

= W jl (C; &2 - &3 - &4 - i2 + r > ^2 + &3 + & 4 + j 2 - -, fa + ir, fci - ir) 



x PFj 2 (p; k 3 -ki~k 2 -ji + ^,k 1 + k 2 + k 3 + j 1 -^,k 4 + ir, fc 4 - ir) 



where 



_ 1 r(fci + k 2 + fc 3 + fc 4 + ji + j 2 - I ± ir)r(fc 2 + fc 3 - I ± «j) 
M ~2vr r(fci + A; 2 + k 3 + j x - \ ± ip)T(/c 2 + & 3 + fa + J2 - | ± »C) 

T(fa + ip ± iu;)r(A:i — ip ± iw)T(/c4 + i£ ± iu;)T(/c4 — iC ± ^>) 
X r(±2tw) ' 

Theorem 17.81 can be obtained from Theorem 17.51 bv setting a = i(k 3 + fa + j 2 — |) and using 
Lemma 15.31 Note that the Racah coefficient for the tensor product 7^ (8 tt^ ® 71"^ is needed 
to obtain Theorem 17.81 It can be checked, using the algebra involution 1?, that these Racah 
coefficients are the same as the Racah coefficients for the tensor product 7r^ (g> 7r^ <8> vr^ , which 
are given in Theorem 17. 11 

The Biedenharn-Elliott identities for the tensor products 7r + (8>7r + (£nr + ®7r p , 7r + (g)7r + ®7r- P (g)7r + , 
and 7r + ® 7r + (g) 7r p (8) 7r~ lead to identities equivalent to those in Theorems 17.51 and 17.61 Finally 
let us remark that Theorems 17. 51 1731 17.81 can all be considered as connection coefficient formulas 
between two orthogonal bases in certain Hilbert spaces. 

8. Racah coefficients for U q (su(l, 1)) 

Throughout this section we assume that < q < 1. Racah coefficients for the quantized 
universal enveloping algebra U q (su(l, 1)) can be calculated in the same way as for su(l, 1). We 
use notations which make it easy for the reader to compare with the Lie algebra case. 

8.1. The quantized universal enveloping algebra U q (su(l, 1)). The quantized universal 
enveloping algebra U q = U q (su(l, 1)) is the unital, associative, complex algebra generated by K, 
K^ 1 , E, and F, subject to the relations 

K 2 - K~ 2 

KK~ 1 = 1 = K~ 1 K, KE = qEK, KF = q~ 1 FK, EF - FE = -—. 

q-q 1 

The Casimir element 

" - r T t "K ~ 2 + EF - «"T \ + " 2 + ^ (8-D 
(q 1 - q) 2 (q- 1 - q) 2 

is a central element oiU q . The algebra lA q is a Hopf *-algebra with comultiplication A defined 
on the generators by 

A(K)=K®K, A(E) = K ® E + E ® K' 1 , 

(8 2) 

A(K- 1 )=K- 1 ®K~ 1 , A(F) = K ® F + F ® R- 1 , 

and A can be extended to U q as an algebra homomorphism. Observe here that the coproduct is 
not cocommutative. The *-structure is defined by 

K* = K, E* = —F, F* = -E, (R- 1 )* = K' 1 . 
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The irreducible ^representations of 1A„ have been determined by Burban and Klimyk There 
are, besides the trivial representation, five classes of irreducible ^representations: 

Positive discrete series. The positive discrete series 7r^~ are labelled by k > 0. The represen- 
tation space is £ 2 (Z>o) with orthonormal basis {e n }neZ >0 - The action is given by 

tt+(K) e n = q k+n e n , vr^ 1 ) e n = g -< fc+ »> e n , 



(q- 1 - q)K+{E) e n = q'^^J (1 - g 2 «+ 2 )(l - q ik + 2n ) e n+1 , 

\ (8-3) 

(a' 1 - q)n+(F) e n = -ql'^^l - ^)(1 _ g 4 fe+ 2n- 2) ^ 

(q' 1 ~ q) 2 *£m e n = (q 2 ^ 1 + q l ~ 2k - 2) e n . 

Negative discrete series. The negative discrete series tt^ are labelled by k > 0. The represen- 
tation space is £ 2 (Z>o) with orthonormal basis {e n } n( =z> - The action is given by 

n-(K) e n = g~( fc+ ") e n , TT-iK' 1 ) e n = q k+n e n , 



(q- 1 - q)7T k (E) e n = -q^-n/ {1 _ q2n){1 _ g 4 fc+ 2n-2) 

/ ( 8 ' 4 ) 

{q~ l - q)H(F) e n = q-2- k - n J(l - g 2 «+ 2 )(l - q* k+2n ) e n+1 , 



(q- 1 ~ q?v k m e n = (q 2 ^ 1 + q 1 ' 2 " - 2) e n . 

Principal unitary series. The principal unitary series representations 7r^ e are labelled by 
< p < — and e G [0,1), where (p, e) ^ (0,^). The representation space is £ 2 (Z) with 
orthonormal basis {e n } n€ z- The action is given by 

rf i£ (K) e n = q n+£ e n , ^{K' 1 ) e n = q~^ e n , 

(q- 1 - q)n P p£ (E) e n = q^^ J{\ - q^+2e+2i P+ l ) (1 _ ? 2n + 2 e -2i„+l) 6n+u 

\ (8-5) 

(q- 1 - q)rf e (F) e n = -q^'^ (1 - q 2n+2e+2i P -l ){1 _ g 2n+2 e -«p-l) ^ 

(q- 1 - (?) 2 < £ (0) e n = {q 2i " + g" 2 * - 2) e n . 

For (p, e) = (0, i) the representation splits into a direct sum of a positive and a negative 
discrete series representation: tt p x = 7r^~©7r7- The representation space splits into two invariant 

U '2 2 2 

subspaces: {e n \ n < 0} © {e n | n > 0}. 

Complementary series. The complementary series representations ir^ e are labelled by A and 
e, where e G [0, ^) and A € (— e), or e G (^,1) and A G (— |,e — 1). The representation 
space is £ 2 (Z,) with orthonormal basis {e n } n ^z- The action of the generators is the same as for 
the principal unitary series, with ip replaced by A + ^. 

The fifth class consists of the strange series representations. The actions of the generators 
in the strange series are the same as in principal unitary series with ip replaced by — + c, 
a > 0. We do not need the strange series representation in this paper. 



8.2. Clebsch-Gordan coefficients. The Clebsch-Gordan coefficients for U q that we need in 
this paper are ^-analogues of the (continuous) dual Hahn polynomials, which we needed to 
describe the Clebsch-Gordan coefficients for su(l, 1). 
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The g-Racah polynomial [3] are defined by 

R n (x;a,P,-f,d;q) = 4 if 3 [ ;q,q 

V aq,/jdq,<yq 

where one of the lower parameters is equal to g _Ar for a non-negative integer iV. These are 
polynomials of degree n in in q~ x + j5q x+1 . 

The dual g-Hahn polynomials are g-Racah polynomials with aq = q~ N and f3 = 0. So the 
dual g-Hahn polynomial are defined by 

'_ iy , (8.6) 

where AMs a non-negative integer, and n € {0, 1, . . . , N}. Under certain condition on 7 and 
5 these polynomials are orthogonal on the finite set {0,1,..., N} with respect to a positive 
measure. The orthonormal dual g-Hahn polynomials are defined by 



t n (x) = t n (x; 7, S, N; q) = (jSq^J ^-//^ ^ <*> N > «)> 

and the orthogonality relations read, for < 7 < g _1 and < 5 < g -1 or for 7 > q~ N and 

AT 

^t m (x;7,<5, N;q)t n (x; 7,(5, N;q)w(x; -y,5, N;q) = 8 mn , 
x=0 

^;i,S,N; q ) = (-ir (-w)"-.,"-*"-' (1 " ^ O t — • 

{Q,dq;q)x{7Sq x + 1 ;q) N +l 

The dual g-Hahn polynomials satisfy the following three-term recurrence relation: 
{q~ x + l5q x+1 ) t n (x) = a n t n+ i{x) + b n t n (x) + a n _it n _i(x), 

where 

a n = ^7^(1 - 9 n+1 )(l - 79 n+1 )(l " q n ~ N /5)(l - q n ~ N ), 
b n = q n ~ N + iq n+l + jq n - N + 7 <5g n+1 - 7 g 2 "- iv (l + g). 

The Askey- Wilson polynomials [3], see also |20| . are polynomials in x + x _1 defined by 

_ f q~ n ,abcdq n ~ 1 ,ax,a/x \ 

p n (x) = p n (x;a,b,c,d\q) = a {ab,ac,ad;q) n 4^3 ;g,g . (8.7) 



a&, ac, ad 

By Sear's 4993 transformation formula [111 (III. 15)] the polynomials p n are symmetric in the 
parameters a, b, c and d. Let the parameters a, b, c, <i be such that non-real parameters occur in 
pairs of complex conjugates, and the pairwise product of two real parameters is smaller than 1. 
Then the Askey- Wilson polynomials satisfy the following orthogonality relations, 

-j— [ p m (x)p n (x)w(x)— + y~] Pm(eq k )p n (eq k ) Res = 5 mn h n . 

4m J T X f-' x=eq k X 

\eq k \>l 

Here T is the counter-clockwise oriented unit circle in the complex plane, e is any of the real 
parameters, and 

(x ±2 ;q)oo 



w(x) 
K 



(ax^ 1 , bx^ 1 , cx^ 1 , dx^ 1 ; q)^ ' 

(abcdq n ~ l ; q) n (abcdq 2n ; q)^ 

(q n+1 , abq n , acq 11 , adq n , bcq n , bdq n , cdq n ; 
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Observe that the Askey- Wilson polynomials can be considered as g-analogues of the Wilson 
polynomials. 

The continuous dual g-Hahn polynomials S n are Askey- Wilson polynomials with d = 0; 

S n (x;a,b,c\q) = p n (x; a, b, c, 0\q) = a~ n (ab,ac;q) nS ip 2 ( '® x > a / x . q q j _ (gg) 

\ ab, ac J 

These polynomials are symmetric in the parameters a, b, c. Let us assume that either a, b, c are 
real and \ab\, \ac\, \bc\ < 1, or one of the parameters is real and the other two parameters are 
complex conjugates with positive real part. In this case the polynomials 

/ \ / 7 I \ S n [x;a,b,c\q) 
Sn{x) = s n {x; a, b, c\q) = = 

y/{q, ab, ac;q) n 

are orthonormal with respect to the measure //(■; a, b, c\q) defined by 

f(x)dji(x;a,b,c\q) = — j f{x)w(x;a,b,c\q) h f(eq k )w k (a;b,c;q), (8.9) 



w(x; a, b, c\q) 



fcez> 

\eq k \>l 

(q,ab, ac, be; q)oo(x ±2 ;q) c 



(ax ±1 ,6x ±1 ,cx ±1 ;g) 00 
where e is any of the real parameters a, b, c. If we assume e = a, we have 

w[x; a, b, c\q) 



Wk{a; b, c; q) = Res 

i=aq k X 

l-a 2 q 2k (g- 2 ,bc;q) 00 (a 2 ,ab,ac;q) k ^ i)* g -§fc(*-i)( 2frc)-*- 



I- a 2 (b/a, c/a; g)oo (q, aq/b, aq/c; q) k 

The polynomials s n form an orthonormal basis for the Hilbert space £(/i) consisting of functions 
/ that satisfy f(x) = f{x~ l ) (/z-almost everywhere), with the inner product associated to the 
measure [i. The polynomials s n satisfy the three-term recurrence relation 

(x + x~ s n (x) (x) + b n s n {x) + (x), (8.10) 

where 



a n = V(i - q n+1 )(l - abq n )(l - acq n )(l - bcq n ), 

b n = a + a' 1 - a' 1 (I - abq n ){l - acq 71 ) - o(l - q n ){l - fec^ 1 ). 

First we consider the tensor product tt^ <£> . This representation can be decomposed into 
irreducible representations as 

oo 

<®<=0< +fc2+ ,- ( 8 - n ) 

3=0 

The Clebsch-Gordan coefficients can be obtained by identifying the action of the Casimir op- 
erator on a suitable subspace of £ 2 (Z>o)® 2 with the Jacobi operator corresponding to the dual 
g 2 -Hahn polynomials, analogously to Proposition 12.11 See also ^Hl, jHZj- |22j . [25] for the inter- 
pretation of the (dual) g-Hahn polynomials as Clebsch-Gordan coefficients. Let us define the 
functions 



4jfci-2 „4fc 2 -2 



q ,p;q 

ip 



These functions are orthonormal on {0, 1, . . . ,p} with respect to the weight 
W(x;k 1 ,k 2 ) = q 4 



4k ix -x(x-i)_ \ L q ){q >q ) 



'l _ q 4k 1 +4k 2 +2x-2^ q 2^ q 4k 2 . q 2^ 
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The intertwiner for Q8.11JI can now be defined by 



x=0 
ni+ri2 



+n 2 -xi 



x=0 



The intertwining property of I' can be checked by writing the Clebsch-Gordan coefficients as 
3^2-functions, and using the following contiguous relations: 



(l-b/q)(l-c/q) 
(1 - d/q){l - e/q)tf 
- n ,b/q,c/ 
d/q,e/q 



3^2 



~ n ,b,c 
d, e 



;q,q 



3^2 



3^2 



-n-1 



, b/q, c/q 



d/q,e/c, 



;q,q) , 



(8.12) 



q n ,b,c 
d, e 



[1 - dq n )(l - eq n ) fq- n ,bq,cq 
(l-d)(l-e) 3iP \ dq,eq 



+ 



[l-q n )(l-deq 



n-l 



■ 3^2 



q n ,bq,cq 



(1 — d)(l — e) \ dq,eq 

These relations follow from |151 (2.4)], respectively |15| (2.3)], and the transformation 



3^2 



q~ n Ac 
d, e 



;q,q 



(q n ,b,c;q) r 
(q,d,e;q) n 



■ 3^2 



-n „1— n 



/d, q x ~ n /e deq n 



^- n jb,q x - n jc 



be 



which is obtained from reversing the summation for the 3(/92-series. 



Next we consider the tensor product representation (8> ir k2 - The action of ir^ <S> 7r fc2 (A(fi)) 
on a suitable subspace of ^ 2 (Z>o)® 2 can be identified with a Jacobi operator for continuous dual 
g 2 -Hahn polynomials. We define for p € [0, — 7r/21ng], 



/^2fc 2 -2fci+l±2ip. q2^ 

(q 2 ,q 4k ^q 2 )- P 



xs n (q 2i P;q 2k z- 2k i- 2 r + \q 2k i+ 2k *-\q 2k i- 2k * +1 \q 2 ), p<0, 

I (g2fci-2fc 2 +l±2ip. g2^ p 



(q 2 ,q^;q 2 ) p 



Xs n (q 2i P-q 2k ^ 2k2 + 2 P +1 ,q 2k ^ 2k ^\q 2k ^ 2 ^ +1 \q 2 ), p>0, 



and similarly for q 2lp in the discrete part of the support of the corresponding orthogonality 
measure. We denote the orthonormality measure for the functions F% by p\ (■; ki, k 2 ). So, 
assuming that k\ and k 2 are such that p? x is absolutely continuous, we have 



f{p)dpP x {p; k u k 2 ) 



2?r 



2 In q 



HP) 



(q 2 ,q^,q 4k2 ,q ±4i P;q 2 ) 



(q2k 1 -2k 2 +l±2ip ^ g2k 2 ~2k 1 +l±2ip ^ q2k 1 +2k 2 ~l±2ip.^ ^2^ 



■ dp 
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As in the Lie algebra case we will define the measure /ii(-; fci, k 2 ) to be the measure Pi(-; ki, k 2 ) 
with the largest possible support. We define the unitary operator I[ by 

TV 

2 In q 







/ 21n9 F™ 1 " 2 (p;A;i,A;2)e ni _ n2 d / 9, 

JO 



where n = min{ni,n2}. With this intertwining operator we find the following decomposition, 

cf. QS1, Q3], 



2 In g 



^ ® 7r fa - J ^pM-k 2 dp- 



This can be proved by checking the actions of K, K~ l , E, F, using the contiguous relations 
1)8.1 2|) . If d/ii has discrete mass points, then we must replace n^ kl _ k2 in the discrete mass points 
by a discrete series representation or a complementary series representation. See |13) Theorem 
2.4] for a precise formulation. 

Remark 8.1. Different from the Lie algebra case, we cannot obtain the tensor product de- 
composition of 7r^T <8> vr^ 2 from the decomposition of n£ ® 7r^ , because the coproduct is not 
cocommutative. When we would try our approach to find the decomposition for ir^ (g) 7r^~ , we 
encounter the following problems. Let Uq PP be the Hopf-algebra U q with opposite comultiplica- 
tion A opp . Then Uq PP is isomorphic to U q -i as a Hopf-algebra, where the isomorphism is the 
interchanging of E and F. So 7r^ 8'7r^ 2 (A(r2)) can be identified with the Jacobi operator for the 
continuous dual g _2 -Hahn polynomials. This Jacobi operator corresponds to an indeterminate 
moment problem, therefore <8> tt^ (A(f2)) is not essentially self-adjoint. We refer to the book 
by Akhiezer for more information about moment problems. 

The action of the Casimir in the tensor product representation ir k ^TTp £ on a suitable subspace 
of £ 2 (Z>o) (g> £ 2 (Z,) can again be identified with a Jacobi operator corresponding to continuous 
dual g 2 -Hahn polynomials. We define 

GP(r;k,p,£) = 



( n n . (g 1+2£±2tP ;g 2 )p (J2iT.J2k+2ip f 2k-2ip l-2fc-2p-2ei 2s , r n 

V (g 1+2fe+2£±2ir ;g 2 ) p 21n9 

( 11" / g #fc (g 1+2£±2 ' P ; g 2 )p / J73 2k+2ip 2k~2ip l-2k-2p-2e\J2\ <Q 

\ {q 2+2j , q 4k+4e ~ 2j ; q 2 ) P 

where Tj = ^ — k — e + j for j £ Z. The functions Gn are orthonormal with respect to a measure 
dp? 2 {-;k,p,E), given by 

//" _ 2 In q \ 
f(T)dp 2 (r;k,p,e)= f{r)w 2 {T\ k, p, e) dr + > f{T)w 2 {Tj;k,p,e), 
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where 



1 (a 2 a Ak a x ~' le± ' li P a ±AiT -a 2 ) 

W2[r;K,p, £) — ^ Sq2k+2ip±2iT^q2k-2ip±2iT^ql-2k-2e±2iT.q2\ ' T fc t U ' _ 21riV' 

i _ „2-4fc-4e+4j /„ 4fc+4e-2 „4fc._2N /„2-4fc-4e _l-2£±2ip. 9\ 

Wn ( T -k n f) = q ^ ,q ,q )co [q — ,q >j 

W2\' 3 ^,P,t) l _ q 2-Ak-Ae ( g 4fc+2 £ -l±2ip ;g 2) oo (^3-4^-2^2^2^ 



We see that the number of discrete mass points of the measure dp? 2 tends to inifinity as p — > oo. 
We denote the measure we obtain in this way, i.e. with an infinite number of discrete mass 
points, by dp 2 . The intertwining operator I 2 is now given by 



2 In q 



I' 2 :£ 2 (Z> )®£ 2 (Z) -> J £ 2 (Z)w 2 {T;k,p,e)dT(B®e 2 (Z>o)w2{T j ;k,p,e) 

3 

TT 

e ni ®e n2 ^ / 21n9 Gl\ +n2 {T-k,p,e)e ni+n2 dT + V G™; +n2 (r i; fc, p, e) e ni 
Jo j 

This leads to the following decomposition of the tensor product Tit ® 7r^ 



7T 

2 In g 

/,■ 'V- ~ 



This is proved by checking the intertwining property of I 2 for the actions of the generators of 
U q , using the contiguous relations (|8.12|) . 

The Clebsch-Gordan coefficients for the tensor products ir£ ® ir^ e and tt^ 10 ~ s 



' a,e 



tained in the same way as for tt~£ ® 7r^ e . Formally the results follow from the substitutions 
p i-> — i(A + |), respectively /? i-> — i(— + a). 

In order to find the decomposition of the tensor product representation nfe®^ we introduce 
the algebra involution defined on the generators by 

•&{K) = K' 1 , &(E) = F, 0(F) = E, tfiK- 1 ) = K, (8.13) 

and d is extended to U q as an algebra mophism. From (|8.3[) and (|8.4[) it follows that ir^(0(X)) = 
tt^(X) for X G U q . Here we identify the representation space for nt with the representation 
space for nZ ■ Furthermore, if we define a unitary operator U : I 2 (Z) — ► £ 2 (Z) by Ue n = 
(-l) n e- n . Then by dHUJ) we have U o ix p p _ e o = vr£ £ o [7. 

A direct verification shows that Ao$ = i?cg)0o A opp , so we have 

< £ ® vr fc o A ^ vr+ ® vr£_ e o A o 

Since i?(f2) = f2 and the Clebsch-Gordan coefficients are eigenfunctions of A(S1), the Clebsch- 
Gordan coefficients for the tensor product 7r^ e <g> ir^ can be obtained from tt^ <g> 7r^ e . Then in 
the same way as for the Lie algebra case, see subsection 13.31 we find that the Clebsch-Gordan 
coefficients for 7r? e (8)7r^" are the functions (— l) n+p G n p {c; k, p, — e). Furthermore, from the tensor 
product decomposition for ir~£ ® 7r^ £ we obtain 

_ 7T 

2 In q 



WILSON FUNCTION TRANSFORMS RELATED TO RACAH COEFFICIENTS 



11 



Here the direct sum is over all j £ Z such that j — e' > 0. 

8.3. Racah coefficients. We give explicit expressions for the Racah coefficients for various ten- 
sor products. The method is the same as in the Lie algebra case, therefore we omit the details. 
In this section we need the g-analogues of the Wilson functions and the Wilson polynomials, the 
Askey- Wilson functions and the Askey- Wilson polynomials. 

The Askey- Wilson function transform is obtained by Koelink and Stokman [22]. The Askey- 
Wilson function transform can be considered as a g-analogue of the the Wilson transform II. 
To stress the similarity between both integral transform, we use notations for the Askey- Wilson 
transform similar to the notations we used for the Wilson transform II. 

Let the parameters a, b, c, d, t satisfy the following conditions 

t>0, < Bt(a), Bt(b) < 1, a = q/d, b = c. 
We define dual parameters by 

ab _ ac ~ ad ~ d 



\b\e 



i arg a 



a a a at 



Note that if a and d are real, we have a = y abcd/q. 
We define the Askey- Wilson function by 

<j)\(x;a,b,c,d\q) = 

{qaXx^/d.q/dX^ab^ac^qa/d^q)^ ,~r~ w , ~ x r, ijw ( 8 - 14 ) 

~ sWi(abcA/q; ax, a/x, aX, bX, cX; q, q/dX). 

(a6cA;g)oo 

This is a slightly different function than the function used by Koelink and Stokman in [231 • To 
identify cp\ with the Askey- Wilson function in [2B], let ip\ be the function defined by |2ol (3.2)], 
then 

4>\(x; a, b, c, d\q) = G(x, \)ipx(x; a, b, c, d\q), 

where 

G(x, A) = (qx^ 1 /d, qX^ 1 /d, ab, ac,bc, qa/d; q)oo. 

The normalization we use here corresponds to the normalization of the Wilson function de- 
fined by (|4.4|) . The reason for choosing this normalization is that the Askey- Wilson function 
4>\{x; a, b, c, d\q) is symmetric in the parameters a, b, c, q/d, while the function ift\(x; a, b, c, d\q) is 
not. This allows us to define the Askey- Wilson function transform that is unitary for a param- 
eter domain that is different from the parameter domain used in |23| . namely the parameters a, 
b, c, d may be non-real. 

By Bailey's transformation [111 (III. 36)] we have 

(qx ±x / d,qX ±l / d,ab,ac;q)oo (ax, a/x, aX, a/A 

4> x (x;a,b,c,d\q) = 4993 

(g/aa; gjoo \ ab,ac,ad 

-4-1 +1 / ~ ~ \ (8.15) 

(ax ,dX ,qb/d,qc/d;q)oo I qx/d, q/dx, qX/d,q/dX 

(ad/q;q)oo 4^3 1 q b/d,qc/ d,q 2 / ad 

From this identity it follows that the Askey- Wilson function satisfies the duality property 

(f>\(x; a, b, c, d\q) = 4> X {X; a, b, c, d\q), 

and that <j>\(x; a, b, c, d\q) is symmetric in a, q/d and in b, c. So, from a = q/d and b = c it follows 
that for x, X € lUT, the Askey- Wilson functions are real-valued. From Bailey's transformation 
jl 11 (III. 36)] with parameters (a,b,c,d,e, f) 1— > (dbcX/ q,ax,a/x,aX,bX,cX) it even follows that 
4>\ is symmetric in a, b, c,q/d, but we do not need this here. 
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Let H be the weight function 

H(x) = H(x; a, b, c, d\q) = - — ±1 '±iL a ^ ±1 v 

(ax 311 , ox^ 1 , cx^ 1 , qx^ 1 / a; q) 00 tt{tx ZCL ;q) 

Here 6{y\q) = (y,q/y,q)oo denotes the Jacobi theta function. Let Z? be the infinite discrete set 
defined by 

v = {t q - k | feeZ, tg~ fc > 1}. 
We define a positive measure /i(-) = h{-; a, b, c, d; t\q) by 



C f dr 

f(x) dh{x) = — / f(x)H(x)— + C T f(x)Res 



x£T> 

where C is the normalizing constant given by 



V 



C = (g; q)ocy/0(at; q)6(bt; q)9(ct; q)9(qt/d; q). 
Explicitly we have, for x = tq~ k € V, 

Res (m) = i 

v= x V V J (q,q,at,a/t,bt,b/t,ct,c/tqt/d,q/dt;q)oo 

(a/t,b/t,c/t,q/dt;q) k ^ ^ 2k 2 \ , 2fc fe2+fc 



x /,'"'"' ' 7 ,7^ 1 - g^/^ (abct z q/d) 
(q/at,q/bt,q/ct,d/t;q) k \ J 

Let Tl(a,b,c,d;t) be the Hilbert space consisting of functions / that satisfy f(x) = fix^ 1 ) 
(/i-almost everywhere), with inner product (•, -)"H defined by 



(f,9)n = J f(x)g(x)dh(x). 

The Askey- Wilson function transform Q : Tt(a, b, c, d; t) — > 7i(a, 6, c, d; t) is defined by 

(Qf)W = (f,<f>x)n- 

The Askey-Wilson function transform is unitary, with inverse = CJ, where is the Askey- 
Wilson function transform with dual parameters. For real parameters a, b, c, d this is Theorem 
1 in [23], for non-real parameters this can be proved using the same method as in |23j . In this 
section we show that the Askey-Wilson function with non-real parameters has an interpretation 
as a Racah coefficient. Then a similar reasoning as in subsection 15.21 can be used to proof uni- 
tarity of Q in this case. 

To determine explicit expressions for the Racah coefficients we need the following g-analogue 
of Lemma 14.61 and Lemma l5.3f iii). 

Lemma 8.2. The Askey-Wilson function satisfies: 

(i) (p x {x;auj,by,b/y,quj/a\q) = 

^ (atu,at/u;q) n 2n ( q~ n ,tx,t/x \ ( q~ n , bty , bt / 'y 

D L^—( 2~^ 1 3V?2 , 3^2 ,, x ,, n 

{q,a z ;q) n \ atuj,at/uj J \ abt\,abt/ A 

_ (a 2 , byx, by/x, bx/y, b/yx, abtX, abt/X)^ 



(bty,bt/y;q) c 



(ii) (paq^(x;a,b,c,d\q) = 

(-d) n q^ n{n ~ l) (qx ±l /d, abq n , acq n ,bcq n ; q)^ p„(/i(x); a, b, c, d\q), n € Z> . 
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Proof. The first statement is a special case of Cor. 3. 4]: Let ct = d in 26, (3.21)], then the 
3(/?2-function on the right hand side of |261 (3.21)] and the 3 (/^-function in the coefficient Gj 
both become 2^1-functions that can be summed by the g-Gauss sum ^| (II. 8)]. This gives the 
following formula 



E(ac,bc;q)j _ 2 , (q 3 ,cx,c/x \ ( q 3 ,cty,ct/y 
{q,ab;q)j \ ac,bc J \ act, abet/ a 

(6ty ±1 , cty ±l , abtx^/a', a't/b; q) c 



(tyx^ 1 ,tx ±x jy, del, abet/a' , ab 2 t/a'; q)^ 
x ^Wi{ab 2 t/ a q; bx, b/x, aby/a', ab/a'y, bt/a'; q, a't/b). 

Then the result follows from substituting 

(a, b, c, a' , t) 1— > (a/uj, au>, t, a/A, b), 

and the definition (|8.14|) of the Askey- Wilson function. 

The second statement follows directly from substituting A = aq n , n S Z>o, in (|8.15|) . and 
from the definition (|8.7|) of the Askey- Wilson polynomial. □ 

There is also a g-analogue of Lemma I5.3lf i^ , 

4>x(x;ap,by,b/y,qp/a\q) = <p p (y;a\,bx,b/x,q\/a\q), 

which follows directly from the definition of the Askey- Wilson function, but we do not need this 
identity here. 

First we consider the tensor product representation tt^ (g) ir^ <8> tt^ . The Racah coefficients 
can be defined as the kernel in the (possibly discrete) integral transform that maps the eigen- 
funtions of the product of Clebsch-Gordan coefficients corresponding to (71^ <8> n~£ 2 ) ® to the 
Clebsch-Gordan coefficients corresponding to 7r^ (£> (71"^ (g> 7r^). We define the weight function 
w(j; k\,k 2 , k 3 ,r) on Z>o as the mass of the measure 

W(j;k 1 ,k 2 ) x Mi(r;/ci + k 2 +j,k 3 ), 

on a line r = constant. Also we define the measure v{p; k\,k 2 , k 3 ,r) to be the measure defined 
by 

P2(p;k 2 ,k 3 ) x p 3 (T;k 1 ,p,k 2 - fc 3 ), 

on the line r = constant. The Racah coefficients for the tensor product tt^ ® 71"^ <g> 71"^ can now 
be determined from 

v™ +r (j; h,k 2 )F^(r; k 1 + k 2 + j, k 3 ) = 



U{j,p;k 1 ,k 2 ,k 3 ,r)F/ m n (p;k 2 ,k 3 )G r n (T;k 1 , p,k 2 - fe 3 ) du(p;k 1 ,k 2 ,h,r), 

or equivalently from 

F m n (P'i k 2, h)G r n (T; ki,p, k 2 - k 3 ) = 

00 (8 16) 

^U{j,p;k 1 ,k 2 ,k 3 ,T)v r ^ +r {j;k 1 ,k 2 )F^ 3 (T;k 1 + k 2 + j,k 3 ) w(j;h, k 2 ,k 3 ,r). 

3=0 

Using the orthogonality relations for the Clebsch-Gordan coefficients, their explicit expressions 
as 3^2-series, and Lemma 18.21 we can express the Racah coefficient as a multiple of an Askey- 
Wilson polynomial. 
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Theorem 8.3. The Racah coefficient Uj(p;ki,k2,k 3 ,r) is the Askey-Wilson polynomial 

N Pj{q 2ip ; g 2fc 2~ 2fc 3+l 5 g2fe 2 +2fc 3 -l ) g2fci+2ir^2fci-2ir |^2^ 

where 

27r(? -2.jfc 1( ^j(j-l) ^2fei+2fe 2 +2fc 3 -l+2i±2iT. q 2^ 



(q ikl ; Q 2 )j (q 2 ,q 2 ,q 4k3 ,q 4kl+4k2+ ^,q ±iir ; q 



2\ 



loo 



9(q 



2fci+2fc2-2fc 3 +l+2j±2ir. 2\ 



y / (g 4fc i+ 4fc 2+2j ) ^l-2fci-2fc2+2fc 3 -2j±2ir. ^2^. 

Remark 8.4. There are not many group theoretic interpretations of the Wilson polynomials. In 
contrast with this, there are many interpretations in the context of quantum groups of the q- 
analogues of the Wilson polynomials, the Askey-Wilson polynomials. For example, Koornwinder 
P§] showed that they occur as spherical functions for SU g (2), see also jHJ and HJ. In the 
Askey-Wilson polynomials occur as matrix elements for representations of U q (su(l, 1)), in |25| . 
[T3] and |Hj they occur as Clebsch-Gordan coefficients. 

Next let us consider the tensor product representation 71^ (g) 7r^ e ® ir^. Her we need the 
measure k%,p, e, k 3 ,r), which is the restriction of the product measure 

/i 2 (cr; h,p,e) x /i 2 (r; /c 3 ,cr, -A;i - e), 

for r = constant. The Racah coefficients for the tensor product 7r^ ®7r / ^ £ ®7r^ can be determined 
from 

(-1)™ G"~ r (C; -e)G;(r; fe l5 C, e - k 3 ) = 

f . , ( 8 - 17 ) 

/ U(a,(;h,P,£,h, t)G™ +v (a; k% , p, e)G m r (r ; /c 3 , a, - fci - e) du(a;ki, p,s,k 3 ,r). 

The inverse for this formula is obtained by substituting k\ «-> /c 3 and e 1— > — e. This also leads 
to the following symmetry property for the Racah coefficients: 

U(a,(;h,p,e, k 3 ,r) = U((, a; k 3 , p, -e, h, r). (8.18) 

Using the orthogonality of continuous dual g 2 -Hahn polynomials and Lemma l8,21 we can express 
the Racah coefficient as a multiple of an Askey-Wilson function. 

Theorem 8.5. For r G [0, — 7r/21ng) the Racah coefficient U(a,C;ki, p,e,k 3 ,r) is the Askey- 
Wilson function 

N <t> q2i< {q 2i °; g 2fe l+ 2 ^, ^34-2^^2^3-2^^2-2^+2^1^^ 

where 

27T g(gl+2fci-2fc 3 +2 £ ±2»r ;g 2 ) 

~~ ( n 2 „4fci „4fc 3 „±4ir n l+2e+2k 1 ±2ia n l-2e+2k 3 ±2i( . „2\ ' 

w > y > y > y > y > y > y 700 

Observe that (|8.18|) is the same as the duality property for the Askey-Wilson funtions. In 
the same way as for the Lie algebra su(l, 1) we can show that the Racah coefficients are the 
kernel in a unitary integral transform. This integral transform is equivalent to the Askey-Wilson 
transform. Let I denote the unitary integral transform 

1 : C(u(a;ki,p,e,k 3 ,T)) -» C(u((; k 3 , p, -e, ki, r), 
(Xf)(C) = J f(^)U(a,C;k 1 ,p,e,k 3 ,r)du(a;ki,p,£,k 3 ,T). 
This integral transform is related to the Askey-Wilson transform Q in base q 2 , with parameters 

(a, b, C, d, t) = (q 2 ^+ 2i P, g 2*3+Mr ) q 2k 3 -2ir^ q 2-2k 1+ 2ip ^ g l-2e-2 kl ^ 

by 

(if)(0 = E (M-^ogoM^/KO, 
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where 



A(a) = g 2 )ooj m = (g l + 2 fc3 -2 £± 2<. q 2~ 

E -- 



OO! 



' iql-2e±2ip^ gl-2fci+2fe 3 -2£±2ir. g2\ 



I fql+2e±2ip^ gl+2k 1 -2k 3 +2e±2iT. q2\ ' 

Remark 8.6. (i) Similar to the Lie algebra case, for the values of r corresponding to discrete 
series representations, the Askey- Wilson function in Theorem 18.51 reduces to an Askey- Wilson 
polynomial, and the corresponding integral transform is the integral transform associated to 
these Askey- Wilson polynomials. 

(ii) Koelink and Stokman [21] found the Askey- Wilson function transform from the interpre- 
tation of the Askey- Wilson function as a spherical function on the quantum group SU 9 (1,1). 
The duality property of the Askey- Wilson functions comes from the interpretation of the Askey- 
Wilson transform as a diffence Fourier transform related to a rank 1 double affine Hecke algebra 
|36| . The interpretation as Racah coefficients given in Theorem 18.51 gives a different explanation 
for the duality property. 

Let us also give the Racah coefficient for the tensor product representation 7r^ <8> Tvp° e . 
Let the weight function w(j; k±, k 2 , p, e, r) be the mass at the points j G Z>o of the measure 

W(j;ki,k 2 ) x fi 3 (T-,ki + k 2 +j,p,e) 

on the line r = constant. Moreover, let v{a; k±, k 2 ,p, £,t) be the restriction of the product 
measure 

p 3 (a;k 2 ,p,e) x fj, 3 (r; ki, a, k 2 +e) 

to a line one which r = constant. Then the Racah coefficient for the tensor product 7r^ <8>7T^ <8>7T^ £ 
is determined from 

G r r Z^ n _ n {o-;k 2 ,p,e)G r n (T;k 1 ,a,k 2 + e) = 

OO 

^U{j, a; h,k 2 ,p, e, r)^" m (j; fei, k 2 )Cf r r Z 3 m -j(r] h + k 2 + j, p, e)w(j; k x , k 2 ,p, e, r), 

3=0 

or from 

v r n ~ m (j; h, h 2 )CT r Zi n _ j (r;k 1 + k 2 + j, p, e) = 

J U (j, cr; ki,k 2 ,p, e, T)G r r Zm -n( a '> k 2, P, e)G^(r; k 1} a, k 2 + e) dv(o\ h,k 2 ,p, e, r). 

Using orthogonality relations for the Clebsch-Gordan coefficients and Lemma l8.2| we can write 
the Racach coefficient as an Askey- Wilson polynomial. 

Theorem 8.7. The Racah coefficient for the tensor product tc^ ® vr^ <g> iTp £ is given by 
Uj(a; h,k 2 ,p, e)=K Pj (q 2ir7 ; q 2k *+ 2i P, ^ka-Mp ^+2^ q 2 kl -2ir ^ 



where 



N =2vr 



( l-2k 1 -2k 2 -2e±2iT 2k 1 +2k 2 +2j+2ip±2ir 2k 1 +2k 2 +2j-2ip±2ir . 2\ 



loo 



(q ikl ; i 2 )j {q 2 , ^i+ 4 fe+ 4 i, qi-^&ip, q ±iir. q 2- 



OO 



, ( n l-2k 1 -2k 2 -2e-2j±2ir. n 2\ . 

q q y ( g 4fc 1 +4fc 2 +2i ;g 2) j 

Recall that in Lie algebra case the tensor product ir + ® ir~ (g> tt + corresponds to the Wilson 
function transform of type I. So it would be interesting to find the Racah coefficients for tl q 
for the tensor product ir + (g) 7r~ ® ir + , since the corresponding integral transform would give 
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a g-analogue of the Wilson function transform of type I. The problem for finding these Racah 
coefficients is that we do not know the decomposition of 7r _ (g) vr + , cf. Remark 18.11 

8.4. Biedenharn-Elliott identities. In the same way as in subsection 17.21 we find integral 
formulas or summation formulas for the Clebsch-Gordan coefficients and the Racah coefficients. 
These do not give essentially new formulas but only new interpretations of known formulas. 
Therefore we only write down two of them, which contain many other known formulas as special 
cases or as limit cases. Both these formulas are obtained from Biedenharn-Elliott identities. To 
simplify notations, we replace q 2 by q, q ki by Jq, for i = 1,2,3,4, and we replace q ±lp by /9 ±1 
(and similarly for a, £ etc.). Moreover, in the formulas we will assume for simplicity < ki < 1, 
ji £ Z> and p, a, r,(,(G T. 

We write out the Biedenharn-Elliott identity for ir^ (g) 7r^ <g> 7r^ <g> ir^ . To find this identity 
we also need the Racah coefficient for the tensor product 7r^ (8) 7r^ ® 7r^ . This Racah coefficient 
can be expressed as a multiple of the g-Racah polynomial, cf. j25j . 

R j3 (h; q 4k2 ~\ q 4k3 ~ 2 , q- 2 *- 2 *- 2 , q ^+^+ 2 h+2h-2. ^ 

The Biedenharn-Elliott identity leads the following convolution identity for Askey- Wilson poly- 
nomials. 

Theorem 8.8. The Askey- Wilson polynomials satisfy the following generalized convolution iden- 
tity 

kt n (Mi l)n(kh q)n(klk 2 2 k 2 gn+n-± ]q)n 
q (fc§ , k 2 k 2 q^; q) h (k 2 k 2 q 2 ^ ; q) jl+h - h 

x R j3 (h; k 2 /q, k 2 /q, q-^-\k 2 k 2 q^- l ]q ) 
x p jl+j2 _ j3 (a) k 2 k 3 q h+ ^ //c 4 , k 2 k 3 k 4 q h ~^ , k%T, k\/r \q) 
x PjAP'i hq^/h, fak 4 q~2 ,k 2 a, k 2 /a \q) 
=Pjiifi] hq^/h, k 3 hq~2 , kik 2 rq j2 , k x k 2 q h jr \q) 
x p j2 {a; k 2 p, k 2 /p, k X T, fci/r \q). 

Remark 8.9. Koelink and Van der Jeugt |25l Thm.4.10] found this identity using the interpreta- 
tion of the Askey- Wilson polynomials as Clebsch-Gordan coefficients for continuous basis vectors 
for positive discrete series representations of U q (su(l, 1)) . If we replace k^ by k 3 k 4 in Theorem 
18.81 and then we put k 3 = 0, we obtain a generalized convolution identity for continuous dual 
g-Hahn polynomials which can also be obtained from writing out (|8,l(i|) and using the orthogo- 
nality relation for the dual g-Hahn polynomials. Theorem 18.81 contains many other generalized 
convolution identities for orthogonal polynomials as special cases, see |25j . 

From the Biedenharn-Elliott identity related to the tensor product tt^ (8) 7r^ 2 ® 7r^ e ® 7rjT we 
obtain the following theorem. 

Theorem 8.10. The Askey- Wilson functions and Askey-Wilson polynomials satisfy 
4> a (p; hk 2 ^q j , k 4 r, fc 4 /r, /fak 2 \q) Pj((; k 2 a, k 2 /a, k%r, k\jr \q) 

p j (uj;k 2 £,k 2 /(;,k 1 p,ki/p\q) (f) T (u; k x p, h(, fc 4 /£ qp/k x \q) 

x 4>a{^;k 2 C,hC,k4/C,^q/k 2 \q) 

K 3 (uj ±2 , k x k 2 £qip ±l , k x k 2 q j p ±l / £\ g)oo 

(k 2 £u ±1 ,k 2 uj ±1 /£,k 4: (u> ±1 ,k4U ±1 /(,k 1 pu ±1 ,kiuj ±1 /p;q) 00 u ' 




h + h 

33 
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Remark 8.11. If a = k^q J2 in Theorem I8.1UI and we use Lemma I8.2f ii). two Askey- Wilson 
functions in Theorem I8.1UI reduce to Askey- Wilson polynomials, and thus we obtain a formula 
with the following structure 

Pjl{QPh(P) = j (pri^PjA^Pj^du. 

This identity can also be obtained from the Biedenharn-Elliott identity for the tensor product 
71 ti ® 7r A^ < ^ )7r fc3 ® 7r fc 4 - ^he Askey- Wilson polynomials in this identity can be further specialized to 
continuous dual q-Hahn polynomials, and this gives an identity that is equivalent to (|8.17[) . Using 
the orthogonality for one of the continuous dual g-Hahn polynomials inside the integral then 
leads to the bilinear generating function for continuous dual g-Hahn polynomials |261 Cor. 3. 4]. 
This is the generating function from which Lemma l8.2f ii) is derived as a special case. 

Also, using the orthogonality relations for the Askey- Wilson polynomial inside the integral 
in Theorem I8.1U1 we find the bilinear generating function for Askey- Wilson polynomials |261 
Thm.3.3]. It is remarkable that in j2B] the generating function is proved using only matrix 
elements and Clebsch-Gordan coefficients for positive discrete series representations oiU q . 
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